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Abstract 

Gathering  insight  into  the  theater  distribution  process  can  be  a  complex  task, 
especially  when  estimating  potential  beddown  solutions.  Coming  up  with  a  low  cost 
feasible  mixture  of  cargo  vehicles  that  will  support  distribution  of  military  personnel  and 
goods  within  theater  is  currently  a  high  priority  for  force  flow  analysts  at 
USTRANSCOM.  In  the  past,  analysts  used  a  trial  and  error  simulation  process  that  was 
iterative  and  time  consuming.  Recent  research  has  produced  the  Improved  Theater 
Distribution  Model  (ITDM),  which  presents  a  less  time  consuming,  more  precise  method 
to  estimate  beddown  allocations. 

Improving  on  this  research,  two  linear  programming  methods  are  developed  and 
added  to  the  ITDM  that  reduce  baseline  beddown  approximations.  Because  daily 
operational  cost  and  initial  beddown  cost  is  included,  this  ultimately  provides  a 
realistically  lower  cost  feasible  solution  when  modeling  theater  distribution.  The 
improved  beddown  solutions  generated  from  post-processing  results  of  the  ITDM  can  be 
used  as  baselines  for  further  distribution  analysis.  Within  the  construct  of  the  model, 
precise  set  notation  is  carried  over  from  the  Improved  Theater  Distribution  Model  and 
slightly  altered  to  reduce  the  generation  of  unnecessary  variables  and  constraints  with 
large-scale  problems. 
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I.  Introduction 


Background 

While  many  objectives  exist  within  the  US  military  to  ensure  the  United  States 
(US)  remains  the  world’s  most  prominent  military  force,  one  important  and  essential 
capability  that  is  core  to  its  success  is  the  distribution  of  troops  and  needed  supplies  from 
the  continental  United  States  (CONUS)  to  overseas  and  deployed  locations.  Theater 
distribution  is  defined  as  the  flow  of  personnel,  equipment,  and  materiel  within  a  given 
theater  as  necessitated  by  the  geographic  combatant  commander  to  support  theater 
missions  (Joint  Chiefs  of  Staff,  2010).  United  States  Transportation  Command 
(USTRANSCOM),  which  is  responsible  for  movements  of  troops  and  supplies,  ensures 
that  these  logistical  needs  are  met.  To  meet  these  requirements,  USTRANSCOM  spends 
a  considerable  amount  of  time  analyzing  data,  assessing  simulations,  and  determining 
appropriate,  feasible  mixtures  of  vehicles  to  employ  so  that  distribution  of  necessary 
supplies  and  troops  is  possible.  The  deployment  process  of  good  and  troops  involves  flow 
from  a  point  of  origin  to  a  point  of  need.  This  point  of  need  is  typically  the  point  at  which 
a  requirement  exists. 

USTRANSCOM  breaks  down  this  journey  of  supplies  and  personnel  from  the 
point  of  origin  to  the  final  point  of  need  into  three  legs.  The  first  leg  involves  movement 
from  a  point  of  origin  to  a  Point  of  Embarkation  (POE).  This  is  usually  from  some 
starting  CONUS  base  to  a  second  en-route  CONUS  base.  This  leg  is  known  as 
intercontinental  movement.  The  second  leg  involves  flow  of  goods  from  a  POE  to  some 
Point  of  Debarkation  (POD),  also  en-route.  The  POD  can  be  thought  of  as  the  midway 
point,  and  this  second  leg  typically  encompasses  movement  from  a  CONUS  location  to  a 
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distinct  theater  of  operations.  This  leg  is  known  as  intertheater  movement.  The  final  leg, 
commonly  referred  to  as  intratheater  movement,  involves  flow  from  a  POD  in  theater  to 
the  final  destination,  or  point  of  need,  which  constitutes  the  point  at  which  the  supplies 
are  needed  (Joint  Chiefs  of  Staff,  2010).  Figure  1  is  included  below  to  illustrate  this 
process  of  intercontinental,  intertheater,  and  intratheater  distribution  of  personnel  and 
supplies. 


(Joint  Chiefs  of  Staff,  2010,  p.  12) 

Figure  1.  Illustration  of  the  Three  Legs  of  Distrihution  Process 


Extensive  research  has  been  conducted  on  all  three  legs  of  the  journey,  but  the  leg 


that  often  poses  the  most  challenge  for  USTRANSCOM  force  flow  analysts  is  the 
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intratheater  journey  from  POD  to  the  final  destination.  Distribution  within  this  area, 
ealled  theater  distribution,  typieally  involves  the  movement  of  supplies  from  an  aerial  or 
sea  port  over  a  relatively  short  distanee  to  a  point  in  a  eombat  zone  or  deployed  loeation. 
Not  only  is  it  essential  to  get  these  supplies  to  their  final  destination,  but  they  must  also 
reaeh  their  destination  in  a  timely  manner. 

Every  grouping  of  supplies  constitutes  a  requirement,  and  every  requirement  is 
accompanied  by  time  windows  within  which  it  may  be  picked  up  and  must  be  dropped 
off  at  its  next  destination.  In  addition,  each  requirement  has  differing  due  dates  for  each 
leg  of  its  journey  to  final  destination.  For  example,  for  a  requirement  to  be  dropped  off  at 
its  POD  there  is  a  time  window  that  has  an  Earliest  Arrival  Date  (EAD)  and  a  Eatest 
Arrival  Date  (EAD).  The  EAD  describes  the  earliest  time  that  delivery  of  a  requirement 
can  occur  at  its  POD  and  the  LAD  describes  the  latest  point  at  which  said  requirement 
can  be  delivered  to  its  POD.  This  creates  a  time  window  within  which  each  requirement 
can  be  delivered  on  its  first  leg  of  the  journey.  There  is  also  a  Required  Delivery  Date 
(RDD)  for  the  second  leg  that  must  be  met  for  the  requirement  to  be  considered  on  time. 
The  RDD  is  the  latest  date  at  which  a  requirement  can  reach  its  final  destination  or  point 
of  need.  On  top  of  this  information,  each  requirement  has  an  associated  weight,  measured 
in  short  tons. 

Under  the  current  system,  USTRANCOM  organizes  all  of  this  information  in 
what  is  called  a  Time  Phased  Force  Deployment  Data  (TPFDD)  file.  The  TPFDD 
contains  all  necessary  information  to  ensure  force  flow  analysts  at  USTRANSCOM  can 
perform  appropriate  studies  and  determine  a  mix  of  vehicles  that  will  ensure  on  time 
delivery  of  all  requirements.  One  more  measure  that  should  be  considered  is  the 
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Commander’s  Required  Delivery  Date  (CRD).  This  date  extends  beyond  the  RDD  and 
allows  requirements  to  be  delivered  in  a  window  between  the  RDD  and  CRD.  It  is  the 
absolute  delivery  day,  and  is  included  so  analysts  may  asses  the  impacts  of  a  late  delivery 
(Joint  Chiefs  of  Staff,  2011a).  These  impacts  will  be  discussed  further  in  Chapters  3  and 
4.  Some  sample  data  that  is  usually  included  in  a  TPFDD  is  shown  in  Table  1  below. 
_ Tablel.  Example  of  Data  Included  in  a  TPFDD _ 


Requirement 

POE 

EAD 

LAD 

POD 

RDD 

Destination 

Total  Short  Tons 

1 

KFFO  (WPAFB) 

5 

8 

OAIX  (baghram) 

10 

GHOS 

200 

2 

KDOV(DoverAFB) 

7 

10 

OAKB  (Kabul) 

12 

BHEL 

300 

As  previously  stated,  USTRANSCOM  takes  particular  interest  in  the  final  leg 
from  POD  to  final  destination.  To  meet  war  fighter  needs.  There  is  a  progression  of 
methods  with  which  USTRANSOM  handles  the  analysis  of  delivering  troops  and 
supplies  over  the  last  leg.  Currently,  analysts  utilize  various  simulation  software  and  tools 
in  order  to  determine  a  feasible  solution  that  meets  both  constraints  on  vehicles  for  the 
theater  of  interest,  and  constraints  on  vehicles  selected  for  simulation.  This  method, 
however,  only  considers  finding  a  feasible  mix  of  vehicles  to  deliver  requirements  in  a 
TPFDD,  and  does  not  prescribe  an  optimal  scheduling  based  on  certain  military 
objectives.  Realistically,  military  leadership  will  have  several  objectives  they  want  to 
meet  such  as  minimizing  operational  cost,  minimizing  the  number  of  vehicles  to  meet 
requirements,  and  minimizing  late  deliveries.  The  method  also  involves  an  extremely 
time  consuming,  iterative  process  in  which  the  operations  plan  (OPLAN)  and  TPFDD  are 
continuously  updated  until  a  feasible  vehicle  schedule  is  developed.  Longhorn  and 
Kovich,  while  working  for  USTRANSCOM  discovered  that  this  process  is  inefficient  and 


4 


may  delay  the  formulation  of  an  operations  plan  (OPLAN)  and  subsequently  the  delivery 
of  essential  supplies  to  the  people  who  need  them  downrange.  As  these  delays  could 
negatively  affect  military  operations  and  overall  mission  success,  clearly,  theater 
distribution  acts  as  a  crucial  point  in  the  delivery  process. 

In  response.  Longhorn  and  Kovich  formulate  a  better  system  for  determining  an 
optimal  mix  of  vehicles  to  meet  TPFDD  requirements.  Their  idea  involves  an  integer 
programming  optimization  model  that  essentially  provides  a  feasible  schedule  of  vehicles 
for  force  flow  analysts  at  USTRANSCOM  to  start  with  before  performing  simulation  of 
distribution.  Since  analysts  would  start  with  this  already  feasible  and  optimal  schedule  of 
deliveries,  the  time  consuming  process  of  iteratively  trying  to  determine  a  feasible 
schedule  would  theoretically  be  eliminated.  This  vehicle  mixture  would  be  used  as  input 
for  simulation  and  should  provide  a  feasible  starting  point  for  analysis.  The  model 
described  in  their  detailed  report  is  better  known  as  the  Theater  Distribution  Model 
(TDM)  and  will  be  referred  to  throughout  this  paper  (Longhorn  and  Kovich).  The  issue 
with  Longhorn  and  Kovich’ s  model  lies  in  the  fact  that  it  presents  far  too  many  decision 
variables  and  constraints  to  be  computationally  efficient.  In  other  words,  the  way  the 
TDM  is  formulated  makes  it  more  complicated  and  sparse  than  it  really  needs  to  be. 

To  reduce  the  size  and  complexity  of  the  TDM,  Micah  Hafich,  an  Operations 
Research  (OR)  studying  at  the  Air  Force  Institute  of  Technology  (AFIT),  proposed  the 
Reduced  Theater  Distribution  Model  (RTDM).  The  RTDM,  while  still  an  integer 
programming  mathematical  model,  reduces  the  size  of  the  TDM  by  introducing  sets  that 
are  not  contained  in  the  TDM  while  maintaining  all  of  the  original  characteristics  of  the 
TDM.  This  improved  formulation  greatly  reduces  the  computational  needs  to  generate  the 
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model,  thus  saving  time  and  money  for  foree  flow  analysts.  The  pure  integer 
programming  aspect  of  this  model,  however,  limits  its  ability  to  allow  more  than  one 
requirement  to  be  scheduled  to  a  single  vehicle.  For  example,  it  would  make  operational 
sense  to  put  two  20  ton  requirements  with  the  same  final  destination  on  one  vehicle  with 
a  payload  of  40  tons,  thus  reducing  the  necessity  of  using  2  or  more  vehicles  to  transport 
the  requirements  when  they  could  have  been  moved  with  1 .  To  remedy  this  issue,  Hafich 
further  proposes  the  Improved  Theater  Distribution  Model  (ITDM).  The  ITDM,  a  Mixed 
Integer  Programming  (MIP)  model,  allows  a  single  vehicle  to  be  assigned  to  deliver  more 
than  1  requirement  if  it  is  in  fact  feasible  for  both  requirements  to  be  placed  on  that 
vehicle.  In  addition,  the  ITDM  addresses  issues  with  modeling  lateness  present  within  the 
RTDM.  For  a  model  who’s  objective  is  finding  a  least  cost,  minimum  lateness  vehicle 
mixture  solution,  it  is  essential  that  lateness  of  requirements  be  modeled  correctly  to 
achieve  a  truly  optimal  or  close  to  optimal  solution.  The  ITDM  takes  care  of  both  of  these 
problems  presented  by  the  RTDM. 

In  addition  to  the  formulation  of  both  the  RTDM  and  the  ITDM,  Hafich  produces 
an  excel  based  Decision  Support  System  (DSS)  which  solves  the  MIP  model  of  the 
ITDM  and  proposes  a  feasible  and  optimal  vehicle  mixture  based  on  an  inputted  TPFDD 
file.  The  DSS  is  a  macro  embedded  Microsoft  Excel  VBA  program  that  uses  the  Lingo 
optimization  software  to  determine  the  optimal  mix;  it  then  outputs  that  mix  into  an 
organized  Excel  file  that  force  flow  analysts  at  USTRANSCOM  can  easily  decipher  and 
use  for  further  simulation.  The  RTDM,  ITDM,  and  accompanying  DSS  tool  fill  the  gap 
that  existed  after  Longhorn  and  Kovich’s  formulation  of  the  TDM,  and  gives  analysts  a 
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definitive  means  with  whieh  to  more  easily  determine  a  feasible  schedule  of  vehicles  to 
meet  all  requirements  specified  in  a  given  TPFDD. 

One  important  take-away  with  the  ITDM  lies  in  the  fact  that  a  solution  output  by 
Hafich’s  DSS  can  be  used  to  determine  an  appropriate  vehicle  beddown  based  on 
upcoming  requirements  within  the  theater  of  operations.  A  vehicle  beddown  involves  the 
placement  of  various  vehicles  of  varying  mode  and  type  at  specific  locations  within 
theater.  It  can  be  reasonably  assumed  that  the  goal  when  determining  an  appropriate 
vehicle  beddown  at  a  POD  would  be  to  place  the  minimum  number  of  vehicles  necessary 
to  meet  delivery  requirements  outlined  in  a  TPFDD,  as  there  is  typically  a  high  cost 
associated  with  moving  vehicles  in  theater  and  maintaining  them  while  deployed  there. 
Vehicle  limitations  exist  as  well.  The  outputs  to  Hafich’s  model  can  currently  be  used  to 
estimate  beddown  needs  at  specific  PODs  within  theater,  therefore  providing  an  efficient, 
effective  tool  for  force  flow  analysis  efforts  at  USTRANSCOM.  The  importance  of  the 
model  and  subsequent  DSS  tool  revolves  around  its  ability  to  produce  a  feasible  vehicle 
mix  quickly  and  consistently  without  the  time  consuming  trial  and  error  methods  used  by 
USTRANSCOM  in  years  past.  In  addition  to  producing  feasible  vehicle  mixtures  when 
considering  distinctive  theater  and  vehicle  constraints,  the  DSS  tool  gives  an  optimal 
solution  with  the  main  objective  of  producing  a  least  cost  on  time  solution.  This  research 
by  Longhorn,  Kovich  and  Hafich  provides  much  needed  support  to  USTRANSCOM 
force  flow  analysts,  and  will  help  supply  war  fighters  with  necessities  when  they  need 
them. 
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Research  Purpose  and  Objectives 

The  purpose  of  this  research  is  to  serve  as  an  improvement  to  the  current  cargo 
and  personnel  mobility  planning  practices  at  USTRANSCOM.  The  capabilities  provided 
by  Longhorn  and  Kovich’s  TDM  model,  and  Hafich’s  RTDM  and  ITDM  models, 
although  great  improvements  upon  the  trial  and  error  techniques  used  before,  do  not 
necessarily  provide  a  realistically  optimal  mix  of  vehicles  from  a  beddown  standpoint. 

The  ITDM  provides  a  least  cost  on  time  solution  with  the  objectives  set  forth  in  the 
model,  but  does  not  necessarily  provide  the  real  world  least  cost  solution,  as  the  task  of 
moving  large  amounts  of  various  types  of  vehicles  into  theater  is  usually  accompanied 
with  a  high  setup  cost.  In  other  words,  the  daily  cost  of  operating  and  maintaining 
specific  types  of  vehicles  are  considered  in  the  objective,  while  cost  associated  with 
placing  vehicles  in  theater  is  not. 

Currently,  the  ITDM  outputs  a  list  of  movements  that  are  mostly  on  time.  A 
shortfall  within  the  model,  however,  is  that  it  attempts  to  move  large  portions  of  certain 
requirements  in  one  or  two  days  instead  of  spreading  delivery  of  these  goods  over  several 
days,  using  the  entire  delivery  time  window.  For  example,  consider  a  100  ton  requirement 
that  must  be  transported  from  its  POD  to  final  destination  with  a  5  day  time  window  to 
deliver  these  goods  before  its  RDD.  Consider  further,  that  the  route  from  POD  to  point  of 
need  only  allows  mode  air  for  delivery,  and  C-130s  with  a  payload  of  10  short  tons  are 
prescribed  by  the  ITDM  to  make  the  delivery.  The  solution  to  the  ITDM  will  tend  to 
move  these  supplies  using  a  large  number  of  C-130s  in  1  Day,  or  over  the  course  of  2 
days.  The  1  day  movement  would  require  10  C-130s  while  a  2  day  movement  would 
require  5  C-130s.  Although  these  allocations  of  vehicles  minimize  per  vehicle  cost  to 
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transport  the  requirement  on  time,  they  do  not  aceount  for  the  initial  cost  and  logistical 
burden  of  the  accompanying  C-130  beddown.  With  this  considered,  transportation  of  the 
requirement  using  2  C-130s  over  a  5  day  period  present  a  more  desirable  solution  since  it 
only  requires  a  beddown  of  2  aircraft  in  theater.  The  overall  objective,  then,  will  be  to 
minimize  the  number  of  vehicles  needed  in  a  theater  of  operations  to  meet  TPFDD 
requirements. 

The  first  objective  of  this  research  is  to  test  the  ITDM  and  determine  whether  the 
results  from  this  model  will  be  adequate  for  approximating  vehicle  beddowns  and 
locations  within  theater  at  a  minimum  setup  cost.  The  current  solutions  do  not  tend  to 
accomplish  this. 

The  second  objective  of  this  research  is  to  make  any  necessary  additions  and 
improvements  to  the  ITDM  formulation  and  DSS.  In  other  words,  the  ITDM  must  be 
improved  to  take  into  account  the  cost  of  a  feasible  vehicle  beddown  for  a  particular 
POD.  This  will  ensure  that  the  objective  of  minimizing  cost  truly  matches  reality,  taking 
into  account  both  vehicle  operation  and  beddown  cost.  It  is  noteworthy  that  cost  may  not 
necessarily  be  measured  in  currency.  As  a  result,  other  costs  such  as  a  weighted  penalty 
may  need  to  be  factored  into  the  model. 

The  third  objective  of  this  research  is  to  test  the  improvements  made  to  the  ITDM 
and  determine  if  these  improvements  account  for  a  beddown  cost.  This  should  minimize 
the  number  of  vehicles  necessary  to  meet  all  requirements  outlined  in  a  given  TPFDD. 

Fourthly,  the  research  will  attempt  to  give  force  flow  analysts  the  ability  to  tailor 
the  model  based  upon  changing  objectives  such  as  minimizing  late  deliveries,  changing 
the  penalty  within  the  model  for  each  late  delivery,  minimizing  beddown  cost,  and 
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minimizing  vehicle  operational  costs.  The  improved  system  should  model  all  of  these 
objectives  simultaneously,  or  allow  the  option  to  model  just  one  or  any  combination  of 
them,  thus  accommodating  for  real  world  objectives  of  decision  makers  at  TRANSCOM. 

Lastly,  this  research  will  strive  to  make  improved  beddown  approximations  at 
specific  POD  locations  from  the  model  results,  with  the  purpose  of  lowering  the  logistical 
and  monetary  burden  that  the  current  solution  to  the  ITDM  prescribes. 

The  culmination  of  these  research  objectives  will  be  in  the  greater  scheme  of 
things  a  better,  lower  cost  method  for  USTRANSCOM  to  approximate  beddown  needs  at 
locations  within  theater  to  support  contingency  operations.  The  research  will  improve  the 
ITDM  to  model  reality  more  closely  and  reduce  the  necessity  of  bedding  down  large 
numbers  of  vehicles  at  PODs.  Force  flow  analysts  should  be  able  to  more  efficiently 
estimate  a  beddown,  which  should  result  in  limited  vehicle  assets  being  available  for  use 
elsewhere.  It  will  free  up  vehicles,  and  offer  the  opportunity  for  these  vehicles  to  be  used 
for  other  military  objectives  to  support  worldwide  military  operations.  The  model 
improvements  should  allow  more  flexibility  to  force  flow  analysts  at  USTRANSCOM, 
allowing  them  to  change  objectives  based  on  preferences.  Through  this  research, 
improved  beddown  approximations  and  model  flexibility  should  improve  efficiency 
when  planning  for  operations  in  theater,  and  provide  an  improved  least  cost,  on  time 
model  of  reality  to  meet  war  fighter  supply  and  resupply  needs. 

Organization 

The  remainder  of  this  thesis  contains  four  additional  chapters.  Chapter  II  provides  a 
literature  review  of  airlift  optimization  modeling,  the  Pickup  and  Delivery  Problem  with 
Time  Windows,  and  other  relevant  models  focused  on  theater  distribution  and  beddown 
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approximating.  It  also  discusses  Integer  Programming  and  Mini/Max  Programming. 
Additionally,  the  proposed  ITDM  is  introduced  and  explained  in  detail.  In  Chapter  III,  the 
methodology  utilized  in  this  research  is  discussed.  In  particular,  two  models,  the  GIBR  and 
the  MPBR  are  introduced.  Chapter  IV  shows  the  implementation  of  the  methodology  and 
demonstrates  improvements  over  the  ITDM.  Chapter  V  offers  concluding  remarks  and 
discusses  how  this  work  might  be  extended  with  further  research. 
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II.  Literature  Review 


This  chapter  reviews  pertinent  and  relevant  literature,  the  purpose  of  whieh  is  to 
give  a  general  understanding  and  baekground  on  the  theater  distribution  problem  and 
attempts  at  modeling  beddown  approximations.  The  information  will  provide  the  reader 
with  a  brief  background  on  the  issue,  but  is  not  an  all-eneompassing  review  of  researeh  in 
this  area.  The  foeus  will  be  on  four  speeific  research  interests  in  the  air  mobility  and 
theater  distribution  field;  a  baekground  on  past  theater  distribution  models,  more 
speeifieally  beddown  planning  efforts,  the  fixed  eharge  assignment  problem,  goal 
programming,  and  an  in  depth  investigation  of  the  Improved  Theater  Distribution  Model 
(ITDM)  as  formulated  by  Hafieh.  The  ITDM  represents  the  most  reeent  effort  to  solve 
theater  distribution  and  beddown  planning  issues,  and  is  the  primary  motivation  for  this 
research.  Thus,  a  detailed  explanation  of  this  model  formulation  is  neeessary. 
Background 

Distribution  planning  is  an  important  part  of  the  overall  joint  operation  planning 
proeess.  It  must  inelude  detailed  analysis  and  evaluation  of  the  distribution  networks  and 
funetions  supporting  the  end-to-end  distribution  proeess,  as  well  as  eneompass  the  full 
range  of  aetivities  neeessary  to  plan  for  national  mobilization,  deployment,  employment, 
sustainment,  and  redeployment  requirements  of  forees  and  materiel  (Joint  Chiefs  of  Staff, 
2010).  This  importanee  drives  foree  flow  analysis  efforts  at  USTRANSCOM.  The  U.S 
military  currently  attempts  to  carry  out  mobility  planning  by  using  numerous  simulation 
tools  sueh  as  the  Generalized  Air  Mobility  Model  (GAMM).  Mekinzie  and  Barnes  (2004) 
deseribe  several  of  these  tools  and  their  importance  in  mobility  seenario  planning. 
Simulation,  however,  only  tends  to  deseribe  seenarios  rather  than  preseribe  a  vehiele  mix 
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to  meet  supply  requirements.  Although  these  tools  help  identify  and  deseribe  limitations 
in  the  theater  distribution,  no  currently  used  tool  prescribes  the  number  of  mobility 
vehicles  or  mixture  of  vehicles  to  meet  theater  movement  requirements  (Longhorn  and 
Kovich,  2012). 

Force  flow  analysis  revolves  around  planning  for  mobility  requirements  outlined 
in  a  TPFDD,  and  updating  those  plans  as  requirements  change.  Although  optimization 
techniques  in  the  past  have  been  geared  towards  routing  of  vehicles,  this  is  not  a  high 
priority  for  analysts  because  of  the  extremely  variable  conditions  in  a  combat 
environment.  For  these  reasons,  the  creation  of  individual  vehicle  routes  and  schedules  is 
neither  necessary  nor  desired  for  force  flow  analysis.  Instead,  analysts  simply  desire  a 
baseline  vehicle  mixture  that  will  successfully  support  distribution  operations  (Hafich, 
2013).  Longhorn  and  Kovich’s  RTDM  and  Hafich’s  subsequent  ITDM  represent  two 
reasonably  successful  attempts  at  optimizing  feasible  baseline  vehicle  mixtures  to 
successfully  support  distribution  operations. 

These  efforts  opened  the  door  for  follow-on  research  in  several  areas.  Vehicle 
beddown  planning  for  PODs  in  theater,  the  primary  focus  of  this  research,  is  one  of  those 
areas.  Accurate  beddown  planning,  as  discussed  in  chapter  I,  is  of  utmost  importance  to 
the  U.S.  military  and  represents  a  key  problem  currently  faced  by  USTRANSCOM.  As  a 
result,  the  drive  for  analysts  is  not  only  to  find  feasible  mixtures  of  vehicles  to  meet 
requirements,  but  also  to  minimize  the  number  of  vehicles  needed  to  do  so,  thus  reducing 
beddown  approximations  at  PODs  of  interest. 

Although  research  in  this  specific  area  is  minimal,  Hafich  along  with  the 
endeavors  of  Longhorn  and  Kovich,  provide  a  basis.  This  chapter  will  provide  an 
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overview  of  past  optimization  efforts  and  other  approaches  to  theater  distribution,  several 
beddown  approximation  techniques,  and  will  explain  in  detail  the  mathematical 
formulation  of  the  ITDM  and  its  beddown  scheme.  It  will  also  outline  the  general  form  of 
a  fixed  charge  assignment  integer  programming  problem,  and  Mini/Max  Programming 
techniques,  the  purpose  of  which  is  to  provide  an  overview  of  approaches  used  in  the 
methodology  portion  of  this  thesis  to  improve  upon  the  beddown  issue. 

Theater  Distribution  Optimization 

Several  attempts  to  optimize  theater  distribution  were  made  in  the  1990s. 
Rappoport,  Levy,  Toussaint,  and  Golden  (1994)  came  up  with  an  airlift  planning  tool  for 
Military  Airlift  Command  (MAC),  predecessor  to  the  US  Air  Force’s  Air  Mobility 
Command  (AMC),  now  headquartered  at  Scott  AFB.  The  single  transportation  Mode 
(air)  model  assigned  different  airlift  aircraft  types  and  shipment  days  to  specific 
requirements.  Once  these  assignments  were  made,  the  results  were  preprocessed  and  then 
processed  utilizing  a  heuristic  routing  and  scheduling  procedure  the  authors  called  the 
Airlift  Planning  Algorithm  (APA).  The  linear  programming  model  minimizes  the  costs  of 
matching  capacity  to  differing  requirements.  Although  their  model  matched  vehicle  types 
to  the  shipments  as  a  form  of  preprocessing,  it  does  not  prescribe  a  specific  number  of 
vehicles  necessary  to  support  distribution  within  the  network. 

Early  optimization  models  also  included  THRUPUT  II,  developed  at  Naval 
Postgraduate  School  and  discussed  by  Rosenthal  et  al.  (1997).  THRUPUT  II,  a  linear 
programming  model,  described  the  entire  distribution  network  in  its  formulation.  The 
model  takes  given  inputs  of  cargo  and  passengers  to  be  moved,  available  airfields. 
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aircraft,  and  routes,  and  minimizes  penalties  for  late  deliveries  subjeet  to  appropriate 
physical  and  policy  constraints.  Priorities  are  assigned  to  requirements  when  necessary. 

Rink,  Rodin,  Sundarapandian,  &  Redfern  (1999)  applied  a  shortest  path  algorithm 
to  the  AMC  aireraft  routing  problem.  The  model  deseribed  utilizes  a  double-sweep 
algorithm  to  find  the  k  -  shortest  paths  between  an  onload  Location  and  offload  location 
provided  in  a  TPFDD.  Unlike  THRUPUT  II,  however,  this  model  does  not  consider 
lateness  and  its  associated  penalties.  In  addition,  the  shortest  path  may  not  be  the  best 
path,  as  unpredictable  eonditions  such  as  weather,  and  enemy  forees  and  other  threats 
may  hinder  sueeess  of  the  delivery.  Lastly,  the  shortest  path  does  not  aceount  for 
outloading  and  unloading  eonstraints  within  the  network;  there  is  no  guarantee  that 
enough  resources  will  be  available  at  certain  airfields  along  the  path. 

In  addition  to  the  shortcomings  described  in  these  models,  another  important 
eonsideration  is  left  out.  Within  a  true  theater  distribution  network,  multiple  modes  of 
transportation  such  as  air,  road  and  rail  are  utilized  to  earry  out  delivery  of  goods  and 
personnel.  All  airlift  models  discussed  thus  far  only  consider  Mode  air.  Other  Modes 
must  be  considered  to  have  a  realistic  model  of  a  transportation  network.  As  sueh, 
beddown  considerations  should  also  include  multiple  modes  of  transportation  as  well. 
Pickup  and  Delivery  Problem  with  Time  Windows 

Typieally  in  theater  distribution,  the  TPFDD  outlines  a  delivery  window  within 
which  a  requirement  can  be  pieked  up  and  must  be  delivered  to  its  destination.  The 
TPFDD  gives  a  time  window  for  both  the  pickup  at  a  POD  and  the  delivery  at  a  final 
destination.  As  a  result.  The  problem  of  theater  distribution  that  the  US  military  and 
USTRANSCOM  faees  can  be  related  to  a  problem  known  as  the  Pickup  and  Delivery 
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Problem  with  Time  Windows  (PDPTW).  Solutions  to  the  PDPTW  yield  optimal  routes 
for  vehicles  in  which  demand  is  met  within  the  appropriate  time  windows  while  meeting 
capacity  and  precedence  constraints  (Dumas,  Desrosiers,  &  Soumis,  1991). 

Dumas  et  al.  (1991)  made  an  early  attempt  at  formulating  the  PDPTW  utilizing  a 
column  generation  and  shortest  path  sub  problem.  This  model  utilizes  a  homogeneous 
fleet  of  vehicles.  Other  endeavors  into  the  PDPTW  that  utilize  a  homogeneous  fleet 
include  a  Reactive  Tabu  Search  method  employed  by  Nanry  &  Barnes  (2000)  and  a  set 
partitioning  formulation  described  by  Baldacci,  Bartolini,  &  Mingozzi  (2011). 

The  homogeneous  nature  of  these  models  makes  them  somewhat  obsolete. 

Models  that  consider  heterogeneous  fleets  of  vehicles  are  far  more  realistic  and  useful, 
and  have  been  researched.  An  exact  algorithm  for  solving  the  PDPTW  with  multiple 
vehicle  types  was  formulated  by  Lu  &  Dessouky  (2004).  The  model,  known  as  the 
Multiple  Vehicle  Pickup  and  Delivery  Problem  (MVPDP),  does  not  necessarily  include 
time  windows.  A  more  robust  model,  developed  by  Xu,  Chen,  Rajagopal,  &  Arunapuram 
(2003),  considers,  in  addition  to  multiple  vehicles,  multiple  time  windows,  compatibility 
constraints,  and  restrictions  on  travel  time.  This  is  known  as  the  Practical  Pickup  and 
Delivery  Problem  (PPDP). 

One  important  consideration  to  note  is  that  the  PDPTW  usually  assumes  that 
vehicles  are  placed  at  centrally  located  bases  from  which  vehicles  begin  their  routes.  A 
beddown,  however,  involves  the  placement  of  vehicles  in  a  theater  of  operations  to 
support  transportation,  and  not  exclusively  the  point  at  which  vehicles  begin  their 
delivery  routes,  which  could  be  anywhere,  including  CONUS  locations.  Specifically,  this 
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idea  of  a  beddown  needs  to  be  assessed  for  PODs  in  theater.  This  goal  takes  high  priority 
for  force  flow  analysts  at  USTRANSTOM. 

Tabu  Search  Techniques  to  Solve  Theater  Distribution  Problems 

Some  of  the  most  recent  approaches  to  theater  distribution  problems  involve  the 
utilization  of  Tabu  Search.  Crino,  Moore,  Barnes,  &  Nanry  (2004)  approached  the 
problem  by  employing  Group  Theoretic  Tabu  Search.  This  method  outlines  specific 
scheduling  and  routing  of  multiple  modes  of  vehicles  using  Tabu  Search.  It  takes  into 
account  delivery  of  goods  within  time  windows.  Similarly,  Burks,  Moore,  Barnes,  &  Bell 
(2010)  describe  an  implementation  of  an  effective  Adaptive  Tabu  Search  (ATS) 
methodology  for  the  Theater  Distribution  Problem  (TDP).  This  methodology  evaluates 
and  provides  a  routing  and  scheduling  of  theater  transportation  assets  at  the  individual 
asset  level  to  ensure  Time  Definite  Delivery  (TDD)  for  all  demands  (Burks,  Moore, 
Barnes,  &  Bell,  2010).  It  solves  both  the  problem  of  depot  location  selection  and  specific 
vehicle  routing  to  support  delivery  needs.  Both  of  these  models  dictate  vehicle  routes  as 
well  as  schedules  at  an  individual  vehicle  level,  while  the  optimization  techniques 
discussed  typically  prescribed  only  one  of  these. 

Vehicle  routing  and  scheduling  provide  little  practical  insight  for  force  flow 
analysis  and  beddown  approximation,  however,  as  the  daily  conditions  in  a  wartime 
environment  are  so  variable.  This  variability  of  conditions  causes  a  day-to-day  scheduling 
approximation  to  be  much  less  useful  and  effective  than  a  generalized  approximation  of 
vehicles  to  be  placed  in  theater  to  support  operations.  Thus,  Tabu  Search  does  not  provide 
a  very  useful  solution  for  determining  baseline  multimodal  vehicle  beddown 
approximations. 
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Past  Beddown  Approximation  Techniques 

Mobility  aircraft,  trucks  and  trains  must  often  times  eompete  with  other  types  of 
vehieles  for  spaee  at  PODs  in  theater.  As  a  result,  many  efforts  have  been  made  to  solve 
the  problem  of  approximating  beddown  needs  for  speeifie  types  of  vehieles  at  strategie 
loeations.  These  attempts  involve  both  mathematieally  based  and  management  based 
approaehes.  Hathaway  et  al.  (1987)  developed  a  method  to  make  bed-down  deeisions 
utilizing  an  integer  linear  programming  model  to  generate  eandidate  basing  (beddown) 
deeisions  for  analysis  and  testing.  Onee  determined,  eandidate  beddown  solutions  are 
simulated  in  FLEETLIFT  for  further  evaluation.  This  model  eaptures  the  dynamie  effeets 
of  the  availability  of  material  handling  equipment  (MHE),  limited  airfield  ramp  spaee, 
variable  distanee  between  network  airfield  loeations,  and  variable  eombat  attrition  and 
planning  faetors  sueh  as  limited  airerew  and  limited  aireraft  loading  eapaeity  (Hathaway, 
1987) . 

Zeisler  et  al.  (2000)  took  a  different  approaeh  by  employing  a  greedy  heuristie  to 
solve  AMCs  intratheater  airlift  seenario  as  a  multiple  knapsaek  problem.  Instead  of 
preseribing  speeifie  vehiele  mixtures  to  meet  given  TPFDD  requirements,  a  generalized 
throughput  assessment  is  given  for  a  predetermined  vehiele  mixture  and  assignment 
seheme.  This  involves  a  trial  and  error  proeess  of  testing  different  beddown  mixtures  to 
maximize  throughput.  This  proeess  of  what-if  analysis  is  time  eonsuming  and  ineffeetive 
for  foree  flow  analysis.  As  a  response  to  the  need  for  a  better  system,  Salmond  et  al. 
(2005)  developed  a  deeision  analysis  method  for  air  mobility  beddown  planning 
seenarios.  Instead  of  making  beddown  deeisions  through  manual  lookup,  trial-and-error, 
and  eorporate  knowledge,  this  researeh  proposes  a  deeision  analysis  tool  that  eompares 
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hard  requirements  versus  existing  capabilities  and  through  gap  analysis,  identifies 
infrastructure  requirement  shortfalls  and  associated  costs  to  satisfy  these  shortfalls 
(Salmond,  et  al.  2005).  This  model,  however,  only  outlines  shortfalls  for  specific 
beddown  decisions.  Additionally,  none  of  these  methods  prescribe  a  general  fleet  mix  to 
support  specific  intratheater  operational  delivery  needs  and  thus,  are  of  little  use  to  force 
flow  analysts. 

A  few  management  based  approaches  have  been  utilized  to  aid  in  beddown 
scenario  planning  as  well.  Koewler  et  al.  (2003)  discusses  improvements  made  to  the 
Capabilities  Based  Logistics  Planner  (CBLP),  a  tool  used  by  Air  Force  Studies  and 
Analysis  Agency  to  quickly  estimate  beddown  plans.  This  is  a  homogeneous  capabilities- 
based  approach  that  assesses  changes  in  airfield  logistics  capacity  as  more  or  less  aircraft 
are  beddown  at  specific  airfields.  Although  a  heuristic  is  developed  to  estimate  the 
parking  capability  of  airfields,  this  model  is  based  upon  very  simple  mathematics. 
Pennington  et  al.  (2006)  developed  a  Microsoft  Access  Based  Cost  Estimation  Tool  for 
Beddown  Analysis  (CETBA).  The  cost-based  model  is  intended  to  provide  the  analyst 
with  the  maximum  amount  of  both  quantitative  and  qualitative  input  for  any  potential 
decision  to  quickly  identify  infrastructure  requirement  shortfalls  and  associated  costs  to 
satisfy  those  shortfalls  (Koewler,  2003).  Similar  to  Salmond’s  Decision  Analysis 
approach,  this  tool  involves  the  assessment  of  infrastructure  shortfalls. 

Although  all  of  these  models  provide  insight  into  the  effectiveness,  capabilities, 
and  costs  of  specific  beddown  scenarios,  they  fail  to  provide  information  about  a  feasible 
mix  of  vehicles  necessary  to  support  theater  distribution  operations.  USTRANSCOM 
requires  beddown  planning  as  a  long  term  decision  based  on  the  specific  requirements 
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given  in  a  TP  FDD.  Therefore  short  term  routing,  seheduling,  and  beddown  capabilities 
and  cost  analysis  does  not  provide  practical  results  to  plan  an  in-theater  beddown  to 
support  force  flow. 

Fixed-Charge  Assignments  Problem  Overview 

Although  fixed  charge  type  approaches  have  rarely  been  applied  specifically  to 
the  theater  distribution  problem,  their  applicability  in  this  area  is  evident.  Winston  et  al. 
(1991)  describes  a  fixed-charge  problem  as  an  integer  programming  formulation  where 
there  is  a  cost  associated  with  performing  some  activity  at  a  non-zero  level,  independent 
of  the  level  of  the  activity.  These  formulations  are  typically  applied  to  production  and 
location  problems.  In  the  production  problem  sense,  if  some  product  is  produced,  a  one¬ 
time  production  setup  cost  is  incurred  no  matter  how  many  items  of  that  product  are 
manufactured.  When  applied  to  location  problems,  a  decision  is  made  on  where  to  locate 
various  facilities  such  as  plants,  warehouses,  or  business  offices,  and  a  fixed  charge  is 
associated  with  building  or  operating  the  facility. 

For  the  purposes  of  this  research,  consider  the  classic  Facility  Location  Problem. 
Given  a  set  Given  a  set  L  of  customer  locations  and  a  set  F  of  candidate  facility  sites,  you 
must  decide  which  sites  to  build  facilities  on  and  assign  coverage  of  customer  demand  to 
these  sites  so  as  to  minimize  cost.  All  customer  demand  di  must  be  satisfied,  and  each 
facility  has  a  demand  capacity  limit  C.  The  total  cost  is  the  sum  of  the  distances  Cij 
between  facility  j  and  its  assigned  customer  i,  plus  a  fixed  charge^  for  building  a  facility 
at  site  j.  This  model  can  be  formulated  as  the  following  integer  linear  program  (SAS 
Institute  Inc.  2010): 
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Decision  Variables: 


Let  yj  =  1  represent  choosing  site  j  to  build  a  facility,  and  0  otherwise. 

Let  Xij  =  1  represent  the  assignment  of  customer  i  to  facility  j,  and  0  otherwise. 


Minimize  X  X  +  X  f jy  j  (1) 

Subject  To 


II 

(2) 

^  yj 

i&L 

Vi  G  L,y/  G  F 

(3) 

Y,d,x,i<Cyi 

ieL 

yj^F 

(4) 

X..  g{0,1} 

Vi  G  L,y/  G  F 

(5) 

y,e{o,i} 

yjGF 

(6) 

Model  1.  Generalized  Fixed-Charge  Location  Problem  (GFCLP) 

The  objective  function  seen  in  (1)  minimizes  cost.  The  constraint  at  (2)  ensures 

that  each  customer  is  assigned  to  exactly  one  site.  Constraint  (3),  known  as  the  linking 

constraint,  forces  a  facility  to  be  built  if  any  customer  has  been  assigned  to  that  facility. 

Lastly,  constraint  (4)  enforces  the  capacity  limit,  C  at  each  site.  The  beddown  problem 

can  be  related  to  the  GFCLP,  only  in  terms  of  general  integers  rather  than  binary.  Some 

of  the  basic  model  concepts  will  be  utilized  in  Chapter  III. 
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Mini/max  Goal  Programming  Overview 


Mini/max  goal  programming  is  typically  used  to  solve  real  world  problems  with 
multiple  and  often  times  eompeting  objeetives.  This  method  attempts  to  meet  some  goal 
or  set  of  goals  rather  than  just  minimize  or  maximize  some  objective,  as  traditional  math 
programming  models  do.  When  integer  variables  are  introdueed  into  a  goal  programming 
model,  it  becomes  an  integer  goal  programming  model;  these  models  ean  contain  zero- 
one  integer  deeision  variables,  general  integer  variables,  or  a  eombination  of  both. 

Several  objeetives  can  be  utilized  in  goal  programming  formulations.  Typieally,  goal 
programming  seeks  to  minimize  the  sum  of  the  deviations  from  all  goals.  Ragsdale  et  al. 
formulates  this  objeetive  as:  Minimize  '^{d^  ’  where  Jr  and  represent  the 

i 

negative  and  positive  deviations  respectively  from  eaeh  goal  i.  One  speeific  formulation 
deals  with  the  Mini-Max  objeetive,  and  is  typieally  formulated  as  seen  in  model  2. 


Minimize:  Q 

(V) 

Subjeet  to 

d-<Q 

(8) 

d:<Q 

(9) 

d-<Q 

(10) 

Ete... 

Model  2.  Mini-Max  Goal  Programming  Formulation 
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Where  Q  represents  the  maximum  deviation  desired  from  eaeh  goal,  i  =  1,2,. .  .and 
di  represents  those  deviations  from  each  goal,  i.  Constraints  (8),  (9),  and  (10)  ensure  that 
no  deviation,  either  positive  or  negative  for  each  goal  i,  exceeds  a  set  value  Q.  Note  that 
both  positive  and  negative  deviations  can  be  modeled,  allowing  more  flexibility  when 
setting  goals.  Mini/Max  Goal  programming’s  relevance  to  beddown  planning  within  the 
ITDM  will  be  discussed  further  in  Chapter  III,  Methodology.  Although  the  exact 
formulation  is  not  used,  some  basic  concepts  are  drawn  from  this  model. 

Improved  Theater  Distribution  Model  (ITDM) 

ITDM  Overview. 

As  previously  discussed,  Hafich  (2013)  improved  upon  the  Longhorn  and  Kovich 
(2012)  theater  distribution  model  formulation,  the  TDM,  by  designing  the  RTDM  to 
greatly  reduce  size  and  complexity,  and  subsequently  the  ITDM  to  model  lateness  of 
deliveries  more  realistically.  Since  the  ITDM  represents  the  most  successful  theater 
distribution  modeling  attempt  to  date,  it  is  of  particular  interest  for  this  research.  The 
ITDM  attempts  to  find  an  optimum  allocation  of  requirements  to  vehicles  in  an  on-time, 
least-cost  manner,  just  as  the  RTDM  does.  The  only  difference  being,  that  the  ITDM  is 
formulated  as  a  mixed-integer  linear  programming  model,  while  the  RTDM  is  a  pure 
integer  programming  model.  This  formulation  is  necessary  as  payloads  for  deliveries  vary 
for  each  vehicle  allocation.  Thus,  assigning  the  same  penalty  for  two  late  deliveries  with 
the  same  Type  of  cargo  and  differing  payload  sizes  does  not  make  practical  sense.  The 
late  delivery  containing  more  short  tons  of  delivery  should  be  assigned  a  higher  penalty. 
For  this  reason,  continuous  decision  variables  that  represent  the  number  of  short  tons 
being  delivered  are  introduced  into  the  model.  This  difference  in  formulation  ensures  that 
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late  deliveries  of  requirements  given  by  a  TPFDD  are  measured  on  a  per  short  ton  seale 
rather  than  per  vehiele. 

With  the  ITDM,  users  must  seleet  whieh  modes  of  transportation  and  vehiele 
types  they  wish  to  enter  into  the  model.  The  individual  Modes  meM  will  typieally 
eontain  all  or  some  elements  of  the  set  {Air,  Road,  Rail}.  Vehiele  types  are  seleeted  by 
the  user  to  form  a  set  of  vehiele  Types  K.  Eaeh  vehiele  Type  k^K  is  a  speeifie  vehiele 
(e.g.  C-17)  of  a  single  Mode  m,  and  has  two  input  parameters  assoeiated  with  it.  The  first 
parameter  is  the  daily  eost  of  utilizing  vehiele  Type  k,  bu.  This  eost  eould  be  linaneial  in 
nature,  but  it  may  also  be  utilized  as  an  arbitrary  eost  in  order  to  analyze  the  impaet 
eertain  poliey  deeisions  have  upon  solutions.  The  seeond  parameter  is  pu  ,  the  average 
payload  (measured  in  short  tons)  of  a  vehicle  of  Type  k  (Hafich,  2013).  These  parameters 
are  essential  to  the  model. 

The  ITDM  draws  in  data  from  the  TPFDD  being  used  for  force  flow  analysis. 
Each  TPFDD  will  list  a  set  of  Requirements  n  ^Wwith  Abeing  the  number  of 
requirements  listed  in  said  TPFDD.  Each  Requirement  n  also  has  an  associated  POD 
i  e  /and  Destination  j  e  J  where  1  and  J  represent  the  set  of  all  PODs  and  destinations, 
respectively.  Every  delivery  of  Requirement  n  also  has  an  associated  delivery  Day  v  on 
which  it  may  be  delivered  to  its  final  destination.  The  set  V  comprises  the  set  of  all 
possible  delivery  days  on  which  Requirement  n  may  be  delivered  to  its  specified 
Destination  j.  Each  movement  Requirement  n,  to  be  delivered  from  POD  i  to  Destination 
J,  has  a  requirement  weight  which  is  measured  in  short  tons.  Within  the  model,  it  is 
assumed  that  all  requirements  are  standard  cargo  requirements.  Passenger  requirements 

and  any  potential  restrictions  on  outsize  or  oversize  cargo  are  ignored  (Hafich,  2013). 
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Some  vehicle  Modes  m  may  not  have  a  direct  path  between  POD  i  and 
Destination  j  supporting  that  vehicle  mode.  Thus,  Mij,  the  set  of  all  Modes  m  with  direct 
paths  between  POD  i  and  Destination  j  is  defined  to  account  for  those  vehicle  modes  with 
no  direct  path  between  a  certain  POD  i  and  Destination  j.  this  set  reduces  the  number  of 
variables  created  by  the  model.  Within  the  ITDM,  Km  represents  the  set  of  all  vehicle 
Types  k  which  are  also  of  Mode  m.  additionally,  since  the  TPFDD  outlines  time  windows 
within  which  Requirements  n  may  be  picked  up  at  POD  i  and  delivered  to  Destination  j, 
not  all  days  v  within  the  set  V  are  necessarily  eligible  delivery  days  for  Requirement  n. 

To  reduce  the  number  of  variables  further,  the  set  Nip  of  Requirements  n  that  are  eligible 
to  deliver  from  POD  i  to  Destination  j  on  Day  v  is  defined.  All  of  these  sets  described  are 
decomposing  sets  within  the  model.  These  decomposing  sets  are  easily  determined  with 
preprocessing  and  are  of  great  value  in  reducing  problem  size  by  eliminating  extraneous 
decision  variable  creation  within  constraints. 

The  delivery  time  windows  and  associated  parameters  for  the  BAD  and  RDD 
require  further  discussion.  All  of  this  information  is  also  given  by  data  in  a  TPFDD.  The 
variable  adn  specifies  the  day  in  which  Requirement  n  arrives  at  its  given  POD.  It  is 
assumed  that  Requirement  n  may  not  be  picked  up  for  delivery  until  the  Day  after  adn. 
Thus,  it  is  not  possible  for  this  requirement  to  be  picked  up  until  Day  adn+1.  Similarly, 
the  variable  specifies  the  Required  Delivery  Date  (RDD),  or  the  day  in  which 
Requirement  n  must  is  desired  to  be  delivered  to  Destination  j.  The  RDD  however  is  not 
an  absolute  deadline  for  Requirement  n.  Thus,  requirements  may  be  delivered  beyond 
their  RDD.  The  parameter  qdn  is  defined  as  the  maximum  allowable  extension  days 
beyond  the  RDD  in  which  Requirement  n  can  be  delivered  to  its  given  final  Destination, 
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j.  Requirements  delivered  outside  the  time  window  ereated  by  parameters  adn+1  and  rdn 
ineur  a  per  short  ton  late  penalty,  g,  whieh  is  user  speeified  depending  on  preferenees. 
Within  the  ITDM,  the  penalty  variable  g  aetually  represents  the  late  penalty  per  short  ton 
per  day  delivered  late. 

The  set  V  mentioned  previously  eomprises  all  Days  v  within  the  time  window 
deseribed  by  the  minimum  of  adn+1  and  the  maximum  of  rdn+  qdn.  No  deliveries  within 
a  given  TPFDD  may  be  made  outside  of  this  minimum-maximum  window,  and  these 
extraneous  deeision  variables  should  not  be  ereated  by  the  model.  This  explains  the 
reasoning  behind  defining  the  set  of  valid  Days  V.  Additionally,  the  ITDM  allows 
aggregation  of  requirements  if  they  fall  within  the  same  delivery  time  window,  doing 
away  with  the  need  for  one  vehiele  to  be  assigned  to  eaeh  single  requirement. 

Most  theater  distribution  models  diseussed  thus  far  eapture  the  limitations  on 
daily  outloading  at  PODs  and  unloading  at  destinations.  The  ITDM  is  no  exeeption.  In 
faet,  variables  are  ereated  to  deseribe  the  maximum  number  of  Mode  m  vehieles  that  ean 
be  outloaded  at  POD  i  on  Day  v,  given  by  Oimv,  and  the  maximum  number  of  Mode  m 
vehieles  that  ean  be  unloaded  at  Destination  j  on  Day  v,  given  by  Ujmv  This  allows  the 
user  to  define  outload  and  unload  restrietions  at  loeations  in  theater  based  on  real  world 
seenarios  and  aetualities,  whieh  ean  be  provided  by  experts  in  the  field.  This  allows 
flexibility  as  POD  eonditions  eertainly  ehange  over  time.  Sinee  some  PODs  and 
destinations  do  not  support  eertain  Modes  m,  Omv  and  ujmv  will  take  on  a  value  of  zero  in 
eertain  oases. 

The  deeision  variables  in  the  ITDM  are  of  two  types.  Variables  Xijmkv  deseribe  the 
number  of  vehieles  of  Mode  m.  Type  k  that  are  required  on  Day  v  to  deliver  any 
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requirements  from  POD  i  to  Destination  j.  Notice  that  this  general  integer  variable  is  not 
tied  to  any  one  Requirement  n.  Thus,  the  vehicle  allocations  dictated  by  decision 
variables  Xijmkv  may  embody  the  movement  of  one,  or  many  different  requirements 
(Hafich,  2013).  To  allow  for  aggregation  of  multiple  requirements  on  a  single  vehicle,  the 
decision  variables  ynijmkv  are  introduced.  They  represent  the  number  of  short  tons  of 
requirement  n  delivered  from  POD  i  to  Destination  j  on  Mode  m.  Type  k  vehicle(s)  on 
Day  V.  these  variable  are  inherently  related,  because  for  every  short  ton  of  Requirement  n 
delivered  from  POD  i  to  Destination  j  on  Mode  m.  Type  k  vehicle(s)  on  Day  v,  some 
vehicle  must  be  assigned  to  make  that  delivery.  The  linking  constraints  described  later 
ensure  this  requirement  is  met. 

To  fully  understand  the  nature  of  the  ITDM  and  its  effectiveness  at  reducing  the 
complexity  of  the  Longhorn  and  Kovich  TDM,  its  Function  Derived  Tuple  Sets  should  be 
discussed.  These  are:  W,  VF,  LF,  VR,  VO,  and  VU.  These  tuple  sets  are  derived  from 
seven  binary  set  defining  functions,  included  in  (1 1)-(17)  below. 

ITDM  Functions  (Hafich,  2013) 


A(jl,  u)  =  I 


_  fl,  if  Movement  n  delivered  on  Day  v  would  be  on  —  time 

O'  otherwise 


(11) 


5(n,  r’)  =  I 


_  fl,  if  Movement  n  delivered  on  Day  v  would  be  late 
0'  otherwise 


(12) 


C  (m,  k)  =  I 


_  fl,  if  vehicle  of  Type  k  is  a  also  a  Mode  m  vehicle 
0'  otherwise 


(13) 
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(14) 


D(n,  i,j) 


1,  if  Movement  n  is  to  be  delivered  from  POD  i  to  Destination  j 

O'  otherwise 


E(i,  m,  v) 


1,  if  3  a  requirement  n  that  may  outload  at  POD  i  onto  a 

Mode  m  vehicle  on  Day  v 
0'  otherwise 


(15) 


II,  if  3  a  requirement  n  that  may  unload  at  Destination  y  onto  a 

Mode  m  vehicle  on  Day  v  (16) 

0'  otherwise 


G{i,j,  v) 


1,  if  3  a  requirement  n,  from  POD  i  to  Destination  j 

s.  t.  adn  +  1  <  r’  <  rdn  +  qdn 
0'  otherwise 


(17) 


The  set  of  tuples  in  VV,  where  W  =  {(/,  j,  m,k,v)\  G(i,  j,  v)  ■  C(m,  k)  =  1} 
eorresponds  to  valid  vehiele  variables  that  may  take  on  a  value.  Thus,  a  vehicle 
variable  is  created  only  when  the  5-tuple  (/,  j,m,  k,v)  corresponds  to  a  theoretically 
possible  vehicle  assignment  (Hafich,  2012). 

The  continuous  decision  variables  associated  with  flows  of  goods  and  personnel, 
ynijmkv,  motivate  the  necessity  for  two  tuple  sets  in  order  to  reduce  the  number  of 
variables.  The  first.  Valid  Flows,  defined  as  VF  =  {{n,  i,  j,  m,k,v)\  A{n,  v)  ■  C(m,  k)  ■  D{n, 
i,  j)  +  B{n,  v)  ■  C{m,  k)  ■  D{n,  i,  j)  corresponds  to  decision  variables  that  are  defined  only 
if  they  are  valid  on-time  or  late  flows.  The  second.  Late  Flows  correspond  to  valid  flow 
decision  variables  that  are  associated  with  late  shipments  in  theater  to  Destination  j. 
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The  three  remaining  tuple  sets  are  Valid  Routes  (VR),  Valid  Outload  (VO),  and 
Valid  Unload  (VU).  Since  only  one  valid  route  exists  for  each  Requirement  n  moving 
from  POD  i  to  Destination  j,  only  a  single  3 -tuple  exists  for  said  requirement  in  the  set 
VR  =  {(n,  i,  j)  I  D(n,  i,  j)  =  1} .  VO  =  {(/,  m,  v)  |  E(i,  m,  v)  =  1}  describes  the  set  of  3- 
tuples  that  are  defined  only  if  Requirement  n  may  outload  at  POD  i  onto  vehicle  Mode 
m  on  Day  v.  finally,  the  set  VU,  defined  mathematically  by  VU  =  {(j,  m,  v)  \  F(j,  m,  v) 
=  1}  is  very  similar  to  the  Function  Derived  Tuple  Set  VO,  with  the  difference  being 
unloading  at  a  Destination  j.  now  that  all  parameters,  sets,  and  decision  variables  have 
been  described,  the  model  formulation  follows.  Table  1  -  Table  4  below  summarize  the 
sets,  parameters,  and  variables  utilized  in  the  ITDM’s  pure  integer  programming 
formulation. 


Table  2.  ITDM  Basic  Sets  (Hafich,  2013) 


Set 

Description 

N 

Set  of  all  Movement  Requirements  n 

I 

Set  of  all  PODs  i 

J 

Set  of  all  Destinations  j 

M 

Set  of  all  vehicle  Modes  m 

K 

Set  of  all  vehicle  Types  k 

V 

Set  of  all  possible  delivery  Days  v 

Mij 

Set  of  all  Modes  m  with  direct  paths  between  POD  i  and  Destination  j 

Set  of  all  vehicle  Types  k  which  are  of  Mode  m 

Niiv 

Set  of  Requirements  n  that  are  eligible  to  deliver  from  POD  i  to  Destination 
j  on  Day  v 
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Table  3,  ITDM  Function  Derived  Tuple  Sets  (Hafich,  2013) 


Set 

Description 

Mathematical  Notation 

VV 

Valid  Vehicle 

{(f  j,  m,k,v)\Gii,  j,v)-Cim,k)  =  l} 

VF 

Valid  Flows 

{(n,  i,  j,  m,k,v)\  A{n,  v)  ■  C{m,  k)  ■  D{n,  i,  j)  +  B{n,  v)  ■  C(m,  k)  ■ 
D{n,i,j)=  1} 

LF 

Late  Flows 

{(n,  i,  j,  m,k,v)\  B{n,  v)  ■  C(m,  k)  ■  D{n,  i,  j)  =  1} 

VR 

Valid  Routes 

II 

VO 

Valid  Outloading 

{(f  m,v)\E(i,m,v)  =  l} 

VU 

Valid  Unloading 

{ij,m,v)\F{j,m,v)  =  l} 

Table  4.  ITDM  Parameters  (Hafich,  2013) 


Parameter 

Description 

h 

Daily  operating  cost  for  a  Type  k  vehicle 

Pk 

Average  payload  of  Type  k  vehicle 

fnij 

Total  weight  (in  short  tons)  of  Requirement  n  to  be  delivered  from  POD  i  to 
Destination  j 

j/j^ 

Day  in  which  Requirement  n  arrives  at  its  given  POD 

rdn 

Required  Delivery  Date  (RDD)  at  the  given  Destination  for  Requirement  n 

qdyi 

Maximum  allowable  extension  days  beyond  RDD  in  which  Requirement  n  can  be 
delivered  to  given  Destination  (with  penalty) 

g 

Late  penalty  per  Short  Ton  late  per  day 

^  imv 

Maximum  number  of  Mode  m  vehicles  that  can  be  outloaded  at  POD  i  on  Day  v 

^jmv 

Maximum  number  of  Mode  m  vehicles  that  can  be  unloaded  at  Destination  j  on 
Dav  V 

^  ijmk 

Daily  cycles  for  a  Mode  m,  Type  k  vehicle  delivering  from  POD  i  to  Destination 

j 

Table  5.  ITDM  Decision  Variables  (Hafich,  2013) 


Variables 

Description 

^ijmk 

Number  of  vehicles  of  Mode  m ,  Type  k  that  are  required  on  Day  v  to 
deliver  requirement(s)  from  POD  i  to  Destination  j 

ynijmkv 

Short  tons  of  Requirement  n  delivered  from  POD  i  to  Destination  j  on  Mode  m , 
Type  k  vehicle(s)  on  Day  v 
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The  mixed  integer  programming  formulation  of  the  Improved  Theater 
Distribution  Model  (ITDM)  follows  below  in  Model  3. 

Minimize  ^  h^ijmkv  +  S  (18) 

(i,j,m,k  ,v)eW  (nj  ,j,m,k  ,v)eLF 


Subject  to 

rd^  +q^ 

3^ nijmkv  ^nij 

M  ij  K,„  v=adfi  +1 

V  (n,i,j)^VR 

(19) 

^ijmk  ^ijmkv  ^imv 

J 

y{i,m,v)  ^VO 

(20) 

^ijmk  ^ijmkv  ^  jmv 

i  Km 

y{j,m,v)^VU 

(21) 

3^ nijmkv  ~  ^ijmkv  ^ijmk  Pk 

% 

y(i,j,m,k,v)GW 

(22) 

3^  nijmkv  —  ^ 

V(n,  i,  j,  m,  k,  v)  e  VF 

(23) 

^ijmkv^{^}^Z^ 

m,  k,v)  &W 

(24) 

Model  3.  Improved  Theater  Distribution  Model  (ITDM)  (Hafich,  2013) 

The  ITDM  represents  a  significant  improvement  upon  previous  theater  distribution 
modeling,  specifically  when  it  comes  to  integer  programming  formulations.  The  objective 
function  at  (18)  utilizes  the  integer  variables  Xijmkv ,  which  describe  the  number  of  vehicles 
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prescribed  by  the  model,  and  the  continuous  flow  variables  ynijmkv,  which  describe  the 
number  of  short  tons  of  a  requirement  delivered  on  a  specific  day.  This  objective  utilizes 
these  variables  to  minimize  not  only  the  cost  of  operating  the  vehicles  prescribed,  but  also 
the  total  penalty  costs  accrued  from  late  deliveries,  measured  in  per  short  ton  late  as 
opposed  to  per  vehicle.  This  per  short  ton  measurement  introduced  by  the  continuous 
variables  certainly  allows  for  more  realistic  modeling  of  theater  distribution  because  it 
does  not  tie  vehicles  to  specific  requirements.  In  other  words,  a  single  vehicle  can  carry 
portions  of  several  requirements  on  one  trip  instead  of  being  tied  to  just  one  requirement. 

Constraints  at  (19)  ensure  that  the  total  sum  of  the  weight  (short  ton  flow 
variables)  associated  with  a  requirement  equals  the  actual  weight  of  that  requirement. 
Constraints  at  (20)  and  (21)  ensure  that  the  user  specified  outloading  and  unloading 
restrictions  at  certain  locations  are  not  exceeded. 

Constraints  at  (22)  link  the  continuous  flow  variables  with  the  integer  vehicle 
variables.  This  constraint  ensure  that  for  all  decisions  corresponding  to  matching 
(ij,m,k,v)  values,  enough  vehicles  are  allocated  to  provide  transportation  capacity  for 
appropriate  requirements  included  as  part  of  those  flows  (Hafich,  2013).  As  mentioned 
before,  this  formulation  allows  vehicles  to  hold  cargo  from  several  requirements  if 
necessary.  In  addition,  it  allows  late  cargo  from  requirements  to  be  delivered  with  on-time 
cargo  from  requirements.  Finally,  constraints  at  (23)  ensure  that  flow  variables  take  on  a 
nonnegative  value,  while  constraints  at  (24)  ensure  vehicle  variables  take  on  a 
nonnegative,  integer  value. 

ITDM  Beddown  Approximation 

As  previously  discussed,  from  the  outputs  of  the  ITDM,  beddown  approximations 
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can  be  made  for  PODs  in  theater.  This  is  aeeomplished  by  post-proeessing  the  solutions. 
Mathematieally,  the  beddown  of  vehieles  of  Mode  m,  Type  k,  needed  at  POD  i  ean  be 
approximated  by:  (Hafieh,  2013) 

(25, 

This  measure  finds  the  maximum  value  of  each  vehicle  Mode  m  and  Type  k 
alloeation  over  all  days  on  the  outputted  delivery  sehedule,  whieh  represents  the  number 
of  this  type  of  vehieles  needed  at  POD  i  to  make  all  deliveries.  This  measure  assumes 
within  the  model  that  for  vehieles  utilized  on  Day  v,  these  same  vehieles  will  be  available 
for  use  on  all  days  following  Day  v.  That  is,  they  complete  a  full  cyele  eneompassing  the 
trip  from  POD  i  to  Destination  j  and  baek  to  POD  i  again.  Vehieles  will  thus  be  ready  for 
use  on  the  day  following  a  delivery.  This  approaeh  will  be  utilized  for  beddown  analysis 
in  the  methodology  and  results  seetions  of  this  researeh. 

ITDM  Conclusion. 

The  ITDM  was  formulated  as  a  response  to  the  TDMs  shortfalls  and  overly 
eomplex  formulation.  Its  purpose  was  to  greatly  reduee  the  size  of  the  model  by 
eliminating  unneeessary  variables  and  eonstraints,  and  to  model  reality  more  effieiently  by 
measuring  flow  of  goods  in  terms  of  short  tons  delivered.  Both  of  these  objectives  were 
met  sueeessfully.  However,  the  ITDM  is  laeking  in  beddown  approximation  effieiency.  As 
stated,  USTRANSCOM  desires  a  tool  to  effeetively  estimate  long  term  beddowns  of  airlift 
aireraft  and  ground  vehieles  at  strategie  theater  loeations.  As  it  stands,  test  runs  for  the 
ITDM  usually  preseribe  a  large  number  of  vehicles  to  aeeomplish  delivery  of  all 
requirements  in  a  given  TPFDD.  This  translates  to  undesirably  large  beddown  estimations. 
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as  costs  associated  with  vehicle  beddown  ean  be  signifieant.  While  the  ITDM  is  eertainly 
an  improvement  on  the  TDM  and  other  beddown  approximation  models  discussed,  there 
is  room  for  further  researeh  in  this  area. 

Conclusion 

Many  of  the  approaehes  to  solving  the  theater  distribution  problem  mentioned  sueh 
as  prior  airlift  optimization  teehniques,  the  PDPTW,  and  Tabu  search,  although  they  have 
expanded  researeh  and  provided  insight,  are  simply  not  what  foree  flow  analysts  at 
USTRANSCOM  eurrently  desire.  To  provide  a  useful  tool  for  approximating  beddown,  a 
model  should  output  vehiele  types  and  the  speeifie  numbers  of  these  vehieles  neeessary  to 
support  operations.  Most  modeling  attempts  prior  to  the  ITDM  foeused  mainly  on  the 
speeifie  routing  and  seheduling  of  vehieles.  As  previously  mentioned,  ehanging  battlefield 
eonditions  hinder  this  approaeh  from  being  of  use  for  analysis  efforts,  espeeially  for 
beddown  analysis.  Thus,  vehieles  and  transportation  eapability  as  output  rather  than  input 
is  mueh  more  useful  for  estimating  beddown  needs.  This  is  preeisely  the  difference 
between  the  ITDM  and  past  modeling  teehniques,  and  the  reason  why  it  is  a  mueh  more 
effeetive  tool. 

Most  of  these  methods  were  formulated  with  the  purpose  of  solving  several  aspeets 
of  the  theater  distribution  problem  simultaneously.  Several  attempts  at  providing  models 
whieh  solve  the  beddown  issue  speeifieally  have  been  made.  Many  of  these  models, 
however,  are  foeused  on  the  effeetiveness,  eapabilities,  and  eosts  of  speeifie  beddown 
seenarios  already  outlined  rather  than  a  mix  of  vehieles  to  support  theater  distribution.  In 
addition,  the  management  based  beddown  approaehes  discussed  focus  mainly  on 
capabilities  and  infrastrueture  shortfalls  given  an  inputted  fleet  of  vehieles  at  speeifie 


34 


bases.  Although  in  some  cases,  this  approach  could  be  helpful,  analyst  at  force  flow 
conferences  are  not  particularly  concerned  with  the  costs  of  infrastructure  shortfalls. 

Although  the  ITDM  provides  exactly  what  USTRANSCOM  needs  in  terms  of 
output,  further  research  into  the  beddown  issue  is  necessary.  As  it  stands,  the  number  of 
vehicles  needed  to  support  theater  distribution  operations  must  be  reduced  in  order  to 
achieve  a  minimal  beddown  approximation.  The  methodologies  outlined  in  this  thesis  aim 
to  improve  the  mixed  integer  programming  formulation  of  the  ITDM  by  providing  a 
reduced  feasible  vehicle  output  that  will  support  TPFDD  requirements.  Details  on  the 
methodology  of  the  improved  model  are  included  in  Chapter  III. 
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III.  Methodology 


Introduction 

This  research  is  organized  into  three  sections  presenting  different  ways  to  modify 
the  ITDM  through  additions  to  the  objective  function  and  the  set  of  constraints  in  order  to 
improve  beddown  approximations.  Firstly,  work  is  done  to  formulate  the  General  Integer 
Beddown  Reduction  (GIBR)  addition  to  the  ITDM.  This  problem  is  centered  on  the 
creation  of  general  integer  variables  within  the  objective  function,  and  necessary 
constraints  are  outlined.  Additionally,  solutions  are  further  improved  by  separately 
applying  the  Mini/Max  Programming  Beddown  Reduction  (MPBR).  This  approach  tends 
to  smooth  vehicle  allocation  solutions  when  the  ITDM  is  applied  to  a  given  TPFDD, 
providing  a  feasible,  reduced  across-the-board  beddown  approximation.  The  MPBR 
formulation  is  tested  individually,  and  results  are  compared  with  those  of  the  original 
ITDM  for  verification  purposes.  Lastly,  MPBR  is  reformulated  with  all  original 
mathematical  notations  to  provide  a  multi-objective  concept  for  force  flow  analysts. 
Finally,  analysis  is  conducted  on  these  proposed  additions  to  the  ITDM. 

Assumptions 

A  large  number  of  the  model  assumptions  are  outlined  by  Longhorn  and  Kovich 
(2012)  and  Hafich  (2013).  Interested  readers  are  referred  to  their  research  for  these 
assumptions.  Several,  however,  are  worth  mentioning  here.  First,  it  is  assumed  that  once  a 
vehicle  is  moved  into  theater  at  a  certain  POD,  it  is  available  for  use  on  all  subsequent 
days.  This  applies  to  all  vehicles  at  all  PODs  within  the  theater  of  interest.  It  is  also 
assumed  that  vehicles  utilized  for  delivery  of  requirements  return  only  to  their  original 
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beddown  location.  That  is,  vehicles  may  not  be  plaeed  at  multiple  PODs,  nor  will  they 
deliver  to  multiple  destinations  within  a  single  trip.  In  reality,  vehicles  may  be  able  to 
deliver  to  multiple  destinations  before  returning  to  a  POD,  but  this  will  not  be  considered 
for  the  purposes  of  this  researeh.  Additionally,  beddowns  at  destinations  are  not 
considered  with  this  researeh.  It  is  assumed  that  the  cost  of  placing  a  vehicle  in  theater  is 
mueh  greater  than  this  vehicle’s  daily  operating  eost,  whieh  generally  holds  true.  Almost 
trivially,  it  is  assumed  that  when  a  vehiele  is  preseribed  to  perform  delivery  of  a 
Requirement  n,  this  vehicle  must  already  be  loeated  at  the  originating  POD.  Thus,  the 
number  of  vehieles  of  a  certain  Mode  m.  Type  k  must  be  sufficient  to  make  all  deliveries 
outlined  in  the  ITDM  solution.  Simply  put,  distribution  of  requirements  is  not  possible 
unless  the  neeessary  vehicles  are  beddown  at  eaeh  POD  eonsidered  in  the  TPFDD. 
General  Integer  Beddown  Reduction  (GIBR) 

GIBR  Motivation. 

Recall  that  the  ITDM  seeks  to  provide  an  optimal  cost  of  delivering  all 
Requirements  n  in  a  TPFDD  by  minimizing  the  combination  of  vehicle  operating  costs 
and  late  delivery  costs.  The  ITDM  does  not,  however,  attempt  to  reduce  the  addition  of 
excess  vehicles  on  days  subsequent  to  Day  v.  Current  solutions  tend  to  preseribe  delivery 
of  requirements  utilizing  a  varied  number  of  vehieles,  whieh  often  times  increase  as  the 
time  window  progresses.  Realistically,  once  vehieles  are  beddown  at  a  particular  POD, 
they  should  be  utilized  as  mueh  as  possible  as  long  as  requirements  are  available.  Instead 
of  adding  more  vehieles  of  the  same  or  even  a  different  Type  k,  these  vehicles  should  be 
utilized  over  and  over.  From  a  beddown  standpoint,  efficiency  means  utilizing  vehicles 
repetitively.  The  ITDM  simply  does  not  reeognize  the  profitability  to  be  gained  by 
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utilizing  the  same  vehicles  repeatedly  and  thus  reducing  the  number  of  vehicles  and 
maintenance  personnel  necessary  to  meet  delivery  needs. 

To  illustrate  this  concept,  consider  the  objective  function  of  the  original  ITDM 
from  Chapter  II  outlined  in  (18). 

Minimize  2  +  Z  siv  -  (18) 

,v)eW  ,v)^LF 

Recall  that  the  first  portion  of  the  objective  function  (18), 

(i,j,m,k  ,v)^VV 

attempts  to  minimize  the  cumulative  operating  cost  of  all  vehicles,  xtjmkvhy  utilizing  the 
daily  vehicle  operating  cost,  bk  ,  for  a  vehicle  Type  k.  Daily  operation  costs  for  a  vehicle 
of  Type  k  are  constant  regardless  of  day.  Thus,  bk  is  the  same  on  Day  v  as  it  is  on  Day 
v+1,  v+2,  and  all  subsequent  days.  To  reduce  beddown  increases  and  to  ensure  that 
vehicles  are  utilized  efficiently  once  in  theater,  some  penalty  should  be  incurred  for 
increasing  vehicles  of  Mode  m.  Type  k  after  Day  v  within  a  TPFDD  time  window.  For 
example,  suppose  v  u-ji-Am-ciso-i  is  prescribed  by  the  ITDM  as  10.  That  is,  10  Mode  AIR, 
type  C-130  are  required  to  deliver  any  available  requirements  from  POD  il  to  Destination 
jl  on  Day  1.  Now,  since  these  10  C-130s  are  flown  from  POD  il  on  day  1,  they  should  be 
utilized  at  POD  il  as  much  as  possible  on  day  2,  day  3,  and  so  forth  as  long  as 
Requirements  n  exist  to  be  delivered.  If  for  instance  x  n-ji-AiR-ciso-i,  which  is  a  Day  2 
allocation,  is  prescribed  as  12,  the  penalty  is  incurred.  If  this  penalty  is  large  enough,  the 
result  should  be  reduction  or  elimination  of  vehicle  beddown  additions.  This  reutilization 
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of  vehicles  should  result  in  more  efficiency  on  Day  v+i  and  subsequent  days,  and 
ultimately  require  a  smaller  beddown. 

Recall  also  that  the  second  part  of  the  objective  function  at  (18), 

(nj  ,J,m,k  ,v)sLF 

seeks  to  minimize  the  costs  incurred  by  delivering  portions  of  requirements  late.  If 
anything,  this  drives  the  model  to  deliver  requirements  quickly  to  avoid  incurring  any  late 
penalties,  g,  especially  if  g  is  large.  The  issue  is  that  the  ITDM  might  avoid  increasing 
beddown  of  vehicle  Mode  m.  Type  k  at  POD  i,  but  nothing  in  the  objective  function 
pushes  it  to  do  so.  Since  neither  portion  of  the  objective  does  this,  solutions  often  show 
increases  in  vehicles  on  days  subsequent  to  Day  v.  The  ITDM  lacks  a  key  portion  in  the 
objective  function  that  will  help  reduce  the  number  of  beddown  increases  once  a  beddown 
is  set.  The  GIBR  addition  to  Hafichs  ITDM  addresses  inefficient  use  of  extra  vehicles  of 
Mode  m,  Type  k  once  a  beddown  is  estimated  on  Day  v.  It  also  ensures  that  quicker 
beddown  approximation  can  be  taken  directly  from  vehicles  prescribed  by  the  ITDM  on 
Day  V,  the  first  day  within  a  given  time  window  {adn+1  to  rdn+  qd„  ),  instead  of 
calculating  them.  The  next  subsection  will  explain  concepts  developed  for  the  GIBR 
before  addition  to  the  ITDM  model  formulation  is  given. 

GIBR  Overview. 

To  accomplish  efficient  utilization  of  vehicles  it  must  be  ensured  that  once  a 
beddown  of  vehicles  is  estimated  for  a  certain  Day  v,  that  for  every  subsequent  Day  (v+1, 
v+2...),  these  same  vehicles  are  utilized  to  deliver  more  Requirements  n  or  portions  of 
requirements.  This  ensures  reuse  of  vehicles  that  are  already  in  theater,  thus  reducing  the 
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necessity  of  adding  more  to  a  given  PODs  beddown.  In  general,  a  beddown  should  be  a 
long  term  solution  that  supplies  enough  vehicles  to  meet  all  delivery  needs  outlined  in  a 
given  TP  FDD. 

The  first  improvement  approach  outlined  in  this  work  is  an  addition  to  the  ITDM 
called  the  General  Integer  Beddown  Reduction  (GIBR).  The  GIBR  involves  the 
introduction  of  decision  variables,  Zimkv,  which  represents  the  number  of  vehicles  of  Mode 
m.  Type  k  that  are  required  to  be  beddown  at  POD  i  on  Day  v  in  order  to  move  any 
eligible  requirements.  These  variables  stem  directly  from  the  ITDM  integer  variables, 

Xijmkv,  and  define  precisely  what  this  research  seeks  to  improve.  As  mentioned,  one  way  to 
help  ensure  that  in-theater  vehicles  are  utilized  on  subsequent  days  is  to  define  a  penalty, 
or  cost  for  situations  when  the  number  of  beddown  vehicles  Zimkv  is  larger  than  the  number 
of  vehicles  Zimk(v+i)  on  a  previous  day.  This  can  be  represented  mathematically  by  positive 
values  of  the  expression  (z.^  •  The  GIBR  has  a  penalty  for  this  situation  defined 

by  Cmk,  which  is  the  cost  of  bedding  down  an  extra  vehicle  of  mode  m,  type  k  at  any  POD 
on  any  day.  It  is  important  to  note  here  that  Cmk  is  constant  for  all  PODs  i.  When  added  to 
the  objective  of  the  ITDM,  the  penalty  for  increasing  a  beddown  can  be  shown  as 
^mki^imkv  ~ ^imk(v-\))  '  Sincc  thc  objectivc  seeks  to  minimize  values,  when  Cmk  is  large 

enough,  the  ITDM  tends  toward  minimizing  the  value  of  .  Thus,  with  this 

extra  portion  to  the  objective  applied  to  the  ITDM,  increases  in  the  beddown  of  vehicles  of 
Mode  m.  Type  k  on  days  subsequent  to  Day  v  are  reduced  and  possibly  eliminated. 

Because  of  its  straightforward  and  elegant  nature,  the  entire  initial  formulation  of 
the  ITDM  is  preserved  and  utilized  in  formulation  of  the  MPBR.  Readers  are  encouraged 
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to  review  Chapter  II  of  this  researeh  for  a  thorough  explanation  of  the  ITDM  and  its 
mathematieal  components.  It  is  important  to  realize  that  the  GIBR  is  not  a  new 
formulation  of  the  mixed  integer  programming  model.  It  is  primarily  an  addition  to  the 
objective  function  and  inclusion  of  three  new  constraints  involving  relationships  between 
the  ITDM  defined  decision  variables  Xijmkv  and  the  GIBR  defined  decision  variables  Zimkv 
All  original  aspects  of  the  ITDM  including  parameters,  decision  variables,  set  defining 
binary  functions,  basic  sets,  and  function  derived  tuple  sets  remain  the  same.  In  fact, 
several  of  these  aspects  are  utilized  when  defining  the  formulation  of  the  GIBR  addition. 
The  three  sets  utilized  within  the  GIBR  are  the  function  derived  tuple  sets  of  Valid 
Vehicles,  W,  and  Valid  Outloading,  VO,  and  the  basic  set  Km.  Recall  that  the  in  order  to 
determine  W,  the  ITDM  describes  the  set  defining  binary  functions  in  (13)  and  (17). 
These  functions  are  included  below. 


C(m,k)  — 


1,  if  vehicle  of  Type  /c  is  a  also  a  Mode  m  vehicle 
O'  otherwise 


(13) 


fl,  if  3  a  requirement  n,  from  POD  i  to  Destination  j 

s.t.  ad^  +  1  <  u  <  +  qd^  (IV) 

0'  otherwise 

These  binary  functions  are  crucial  in  the  creation  of  vehicle  decision  variables 
within  the  ITDM,  which  populate  the  set  W.  The  set  of  tuples  in  W, 
where  W  =  {{i,j,  m,  k,  v)  \G{i,  j,  v)  ■C{m,  k)=\},  includes  those  tuples  which  correspond 
to  valid  vehicle  variables  that  may  take  on  value  within  the  mixed  integer  program.  The 
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set  of  vehicle  decision  variables  V(/,  j,m,k,v)  is  used  to  describe  summations  in  the 
new  objective  function  as  well  as  which  constraints  the  updated  model  includes  as  valid. 
The  set  VO  =  {(/,  m,  v)  |  E{i,  m,  v)  =  1}  describes  the  set  of  3-tuples  that  are  defined  only 
if  requirement  n  may  outload  at  POD  i  onto  vehicle  Mode  m  on  any  day.  Detail  on  how 
this  set  is  derived  is  given  in  (15)  below. 

fl,  if  3  a  requirement  n  that  may  outload  at  POD  i  onto  a 

Mode  m  vehicle  onDayu  (15) 

O'  otherwise 

In  addition,  the  GIBR  employs  one  of  the  ITDM  basic  sets,  Km,  the  set  of  all 
vehicles  that  are  of  Mode  m.  Type  k.  The  set  Km  is  included  for  describing  summations  in 
the  new  constraints  and  the  objective  function.  It  is  important  to  note  the  beddown 
decision  variables  Zimkv  need  not  be  defined  by  specific  POD  i.  They  can  also  be  described 
across  the  entire  theater  of  operations  including  all  PODs  within  a  TPFDD  utilized  for 
analysis.  In  other  words,  if  analysts  are  interested  in  estimating  beddowns  by  Location  i, 
the  variables  Zimkv  would  be  included  in  the  model.  Further,  if  they  are  interested  in  the 
estimation  of  a  theater  wide  beddown,  the  variables  Zmkv  are  included.  This  describes  the 
number  of  vehicles  of  Mode  m.  Type  k  that  are  required  to  be  beddown  in  a  specific 
theater  (all  PODs  in  TPFDD)  on  Day  v  to  deliver  any  eligible  requirements. 

Several  different  variations  to  the  formulation  exist  such  as  including  variables 
Zimv,  which  represents  the  number  of  vehicles  of  Mode  m  required  to  be  beddown  at  POD 
i  on  Day  v  to  deliver  eligible  requirements.  There  is  also  Zim,  which  represents  the  number 
of  vehicles  of  Mode  m  required  to  be  beddown  at  POD  i  on  any  day  to  deliver  any 
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eligible  requirements.  Redefining  these  deeision  variables  depends  on  analysts  desires. 
Note  that  as  the  deeision  variables  are  indexed  differently,  the  penalty  parameter  must  be 
indexed  differently  as  well.  For  the  purpose  of  this  researeh,  only  the  formulation 
ineluding  the  deeision  variables  Zimkv  and  eost  parameter  Cmk  is  ineluded  and  deseribed  in 
this  ehapter.  The  ITDM  remains  a  mixed  integer  programming  model,  and  any  additions 
made  to  it  are  outlined  and  deseribed  in  Models  3  and  4  below. 

General  Integer  Beddown  Reduction  (GIBR)  Addition. 

Initial  ITDM  Formulation. 

The  ITDM,  whieh  is  the  initial  model,  is  shown  below  in  Model  3. 

Minimize  ^  ^  siv  -  (18) 

(i,J,m,k  ,v)  eVV  (n,i ,  j,m,k  ,v)  ^  LF 


Subjeet  to 

rd  +q 

SEE  y  nijmkv  y(n,i,j)^VR  (19) 

M  jj  v=ad„  +1 

^ijmk  ^ijmkv  -^imv  \f{i,m,v)^VO  (20) 

J  K„ 

'^ijmk^ijmkv  -^jmv  (21) 

/  Km 

'^ynijmkv  -^ijmkv%mkPk  V (i,  j ,  M,  k ,v)  ^  W  (22) 

Nil, 

ynijmkv  ^  0  V(n,  i,  j,  m,  k,  v)  e  FF  (23) 

^ijmkv  e  {0}  ^  v(/,  j,  m,  k,v)  (24) 


Model  3.  Improved  Theater  Distribution  Model  (ITDM)  (Hafieh,  2013) 
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ITDM  With  GIBR  Addition  (By  Location  i,  Mode  m,  Type  k,  Day  v). 


The  General  Integer  Beddown  Reduetion  variation  of  the  ITDM,  which  describes 
beddown  decision  variables  by  Location  i  is  formulated  below  in  Model  4.  Tables  6  and  7 
outline  new  parameters  and  decision  variables.  Any  additions  (26)-(29)  to  the  ITDM  are 
followed  by  an  asterisk. 


Table  6.  GIBR  Parameters 


Parameter 

Description 

^mk 

Penalty  for  bedding  down  extra  vehicle  of  Mode  m.  Type  k  at  any  POD  on 
any  day 

Table  7.  GIBR  Decision  Variables 


Variable 

Description 

Zimkv 

Number  of  vehicles  of  Mode  m.  Type  k  that  are  required  to  be  beddown  at 
POD  i  on  Day  v  to  move  any  eligible  requirements 

GIBR  Formulation. 

The  GIBR,  which  is  the  initial  ITDM  model  with  additions,  is  shown  below  in 
Model  4. 
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Minimize  2  +  S  siv  -  + 

(i,J,m,k  ,v)  ^W  (n,i,j,m,k,v)^  LF 

^  c  Tz  — z  ^  ^ 

mk\'^imkv  ^imk{v—\)^ 


{i,m,v)€VO 


Subject  to 

rd^  +q^ 

EZE  y  nijmkv 
M,i  K,„  v=ad„  +1 


=  r  ■■ 

mj 


^ijmk  ^ijmkv  — 


J  K„ 


^ijmk  ^ijmkv  ~  ^  jmv 


V  (n,i,j)^VR 


\f{i,m,v)  gVO 


y{j,m,v)^VU 


(26) 


(19) 

(20) 

(21) 


/  /r„ 


3^ nijmkv  ^ijmkv  ^ijmk  Pk 


'‘ijv 


y  jc  <7 

ijmkv 

j\{i,j  ,m,k,v)eW 


imkv 


y{i,j,m,k,v)^W 


V(/,  m,  v)  e  VO,  (Vk  e  ^„,) 


* 


(22) 

(27) 


z  ^  z 

^imkv  ^imk{v—\) 


y  nijmkv  —  ^ 

e{0}uZ" 


y(i,m,v)GVO,  (VkGKJ 


V(n,  /,  j,  m,  k,  v)  e  VF 


V(i,m,v)  gVO,  (Vk  ^Km) 


V{i,j,m,  k,v)  e  W 


* 


(28) 

(23) 
(29) 

(24) 


Model  4.  ITDM  With  GIBR  Addition  (By  Location  i,  Mode  m,  Type  k,  Day  v). 
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The  GIBR  presents  a  significant  improvement  over  ITDM  beddown 
approximations.  First,  addition  of  the  cost  parameter  Cmk  ensures  that  some  penalty  is 
assessed  in  the  objective  function  (26)  when  a  Day  v  beddown,  zimkv,  is  increased  on  any 
day  subsequent  to  v.  If  Cmk  is  set  large  enough,  this  penalty  ensures  that  a  beddown  of 
vehicles  of  Mode  m.  Type  k  will  not  increase  throughout  a  time  window.  Because  a 
beddown  will  not  increase,  the  model  will  utilize  the  vehicles  already  in  theater  more  often, 
providing  improved  vehicle  utilization  and  a  more  consistent  vehicle  beddown.  It  also 
provides  simpler  beddown  estimation,  as  the  maximum  of  vehicles  would  be  the  number  of 
each  type  of  vehicle  prescribed  on  the  first  day  in  the  time  window. 

The  objective  function  is  now  threefold.  It  attempts  to  simultaneously  minimize 
vehicle  usage  costs,  late  deliveries  penalties,  and  beddown  increase  penalties.  Later  this 
research  shows  how  manipulating  the  costs  bt  and  penalties  g  and  Cmk  can  cause  the  model 
to  achieve  certain  desired  solutions.  Again,  this  formulation  of  the  GIBR  is  indexed  by 
POD  i,  rather  than  across  the  whole  theater  of  operations.  Thus,  it  seeks  to  improve 
beddowns  by  specific  location  rather  than  by  theater. 

Three  additional  constraints  are  added  to  the  model  as  well.  Constraints  at  (27) 
ensure  that  whatever  the  number  of  vehicles  prescribed  by  the  model,  the  beddown  is 
sufficient  to  meet  these  needs.  Thus,  beddown  of  vehicles  at  every  POD  i  must  be 
sufficient  to  meet  delivery  needs. 

Constraints  at  (28)  ensure  that  on  Day  v+1  that  there  are  at  least  Zimkv  vehicles 
available.  Any  increase  in  Zimkv  on  days  after  v+1  will  be  penalized  by  the  objective 
function.  That  is,  beddown  of  vehicle  Mode  m  Type  k  will  be  consistent  at  a  Location  i 
throughout  a  given  time  window.  Constraints  at  (29)  ensure  that  beddown  variables  are 
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nonnegative  and  integer.  This  is  neeessary  because  a  partial  vehicle  does  not  make  sense. 
Thus,  a  beddown  must  be  a  positive  integer. 

GIBR  Addition  Conclusion. 

The  GIBR  is  the  first  of  two  main  contributions  to  this  thesis.  The  formulation, 
with  the  use  of  new  parameters  and  decision  variables,  forces  the  ITDM  to  more 
efficiently  utilize  vehicles  over  a  time  window.  The  model  now  attempts  to  minimize 
spikes  in  the  number  of  a  certain  type  of  vehicles  later  in  a  TPFDD  time  window,  and  to 
utilize  vehicles  already  in  theater  to  move  requirements  in  that  TPFDD.  It  also  tells 
analysts  when  the  best  time  to  add  vehicles  to  a  beddown  if  in  fact  it  becomes  necessary 
to  add  vehicles.  Thus,  the  GIBR  provides  beddown  reductions  for  force  flow  analysis  and 
can  help  save  valuable  vehicle  resources  which  may  be  utilized  elsewhere. 

Mini/Max  Programming  Beddown  Reduction  (MPBR)  Addition 

MPBR  Motivation. 

The  Mini/Max  formulation  is  mathematically  equivalent  to  the  GIBR.  Further, 
given  the  same  penalty,  it  will  have  the  same  effect  on  the  ITDM  objective  when 
solutions  are  found.  The  exception  to  this  is  when  there  are  multiple  optimal  solutions.  In 
this  case,  the  exact  same  solution  may  not  be  output.  The  difference  between  the  two 
models,  however,  is  the  number  of  variables  needed  to  formulate  the  linear  program  since 
the  MPBR  problem  is  not  indexed  by  v. 

For  this  reason,  some  addition  that  minimizes  the  maximum  of  the  number  of  vehicles 
with  fewer  variables  is  worth  investigating.  Traditionally,  Mini/Max  Programming  is 
applied  to  problems  such  as  this.  Thus,  an  improvement  known  as  the  MPBR  is  developed 
to  reduce  beddown  approximations  and  vehicle  allocations.  This  addition  is  similar  to  the 
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GIBR  in  that  it  involves  a  single  addition  to  the  ITDM  objeetive  funetion  as  well  as  a  few 
eonstraints. 

Reeall  that  the  ITDM  seeks  to  minimize  the  eost  of  delivery  of  all  Requirements 
n^N  in  a  TPFDD  by  minimizing  the  addition  of  vehiele  operating  eosts  and  late  delivery 
eosts.  The  ITDM  does  not,  however,  seek  to  reduee  the  size  of  the  fleet  neeessary  to 
aeeomplish  this  distribution.  Current  solutions  tend  to  preseribe  delivery  of  requirements 
quiekly  using  a  large  fleet,  instead  of  utilizing  the  full  allowable  time  window  and 
fulfilling  requirements  with  a  smaller  fleet.  The  ITDM  simply  does  not  reeognize  the 
profitability  to  be  gained  by  redueing  the  overall  number  of  vehieles  necessary  to  meet 
delivery  needs.  The  result  is  solutions  and  subsequent  beddown  approximations  that  are 
often  much  larger  than  they  need  to  be. 

To  illustrate  this  concept,  consider  again  the  objective  function  of  the  ITDM  from 
Chapter  II  outlined  in  (18). 

Minimize  X  +  X  s{v (18) 

,v)eW  (n,i,J,m,k  ,v)^LF 

As  previously  discussed,  the  first  portion  of  the  objective  function  (18), 

^ijmkv 

,v)eW 

seeks  to  minimize  the  overall  cost  of  utilizing  all  vehicles  Xijmkv  by  utilizing  the  daily 
vehicle  operating  cost  bt  for  a  vehicle  k.  Daily  operation  costs  for  a  vehicle  Type  k  are 
constant  regardless  of  day.  Thus,  bt  is  the  same  on  Day  v  as  it  is  on  Day  v+1  and  all 
subsequent  days. 
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Consider  a  delivery  sehedule  that  requires  a  Mode  m,  Type  k  vehiele  10  times  with 
an  allowable  time  window  {adn+1  to  rJ„+  qdn  )  of  5  days  for  delivery.  Now,  since  bu  is 
constant  across  all  Days  v,  the  timeline  of  this  delivery  is  irrelevant  as  long  as  all  short 
tons  of  Requirements  ne  A^are  delivered  within  the  allowable  time  window.  Assuming  no 
late  deliveries,  there  is  no  difference  in  operating  cost  for  these  requirements  when 
delivering  over  2  days  using  5  vehicles,  as  opposed  to  delivering  it  over  5  days  using  2 
vehicles.  As  a  result,  the  first  part  of  the  objective  does  not  seek  to  minimize  vehicles  used 
overall,  but  seeks  to  minimize  total  cost  of  trips  made  by  vehicles.  The  true  difference  is 
that  realistically,  a  5  vehicle  beddown  costs  more  than  a  2  vehicle  beddown  of  the  same 
mode  and  type. 

Recall  again  that  the  second  part  of  the  objective  function  (18), 

(n,i  ,J,m,k  ,v)^LF 

seeks  to  minimize  the  costs  incurred  by  delivering  portions  of  requirements  late.  If 
anything,  this  drives  the  model  to  deliver  requirements  as  quickly  as  possible  to  avoid 
incurring  any  late  penalties  g,  especially  if  g  is  large.  The  issue  is  that  the  ITDM  might 
spread  delivery  over  a  given  time  window,  but  still,  nothing  in  the  objective  function 
pushes  it  to  do  so.  The  MPBR  addition  to  Hafich’s  ITDM  addresses  high  cost  beddown 
solutions  with  large  numbers  of  vehicles.  Keep  in  mind  the  formulation  of  the  MPBR 
produces  the  same  beddown  solution  as  the  GIBR,  but  eliminating  subscript  v  in  decision 
variables  of  the  MPBR  decreases  the  number  of  variables  produced  when  building  the 
model.  In  other  words,  these  models  are  simply  two  formulations  that  accomplish  the  same 
solution.  The  next  subsection  explains  concepts  developed  for  the  MPBR  before  addition 
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to  the  ITDM  model  formulation  is  given  and  deseribed. 

MPBR  Overview. 

The  MPBR  introduees  a  new  eost  parameter  and  a  new  set  of  deeision  variables  in 
order  to  smooth  vehiele  solutions  to  the  ITDM.  Instead  of  assessing  vehiele  eosts  based 
solely  on  daily  operation,  a  beddown  eost,  dm ,  indexed  by  vehiele  Mode  m  is  utilized.  As 
deseribed  previously,  these  extra  eosts  provide  a  more  realistic  assessment  of  overall 
costs.  In  addition,  decision  variables  Qm,  which  represent  the  upper  bound  on  the  number 
of  Mode  m  vehicles  to  be  beddown  at  POD  i  on  any  day  within  the  allotted  time  window 
are  introduced.  As  with  the  GIBR,  the  MPBR  is  not  a  new  formulation  of  the  ITDM.  It 
simply  involves  the  addition  of  a  portion  to  the  objective  function  along  with  a  few 
necessary  constraints. 

All  original  aspects  of  the  ITDM  including  parameters,  decision  variables,  set 
defining  binary  functions,  basic  sets,  and  function  derived  tuple  sets  remain  the  same.  In 
fact,  three  of  these  aspects  are  utilized  when  defining  the  formulation  of  the  MPBR 
addition.  The  tuple  sets  utilized  by  the  MPBR  are  the  set  of  Valid  Outloading,  VO,  and  the 
set  of  Valid  Vehicles,  W.  Detail  on  how  these  sets  are  derived  is  given  in  (15)  previously. 
The  basic  set  Km  is  used  in  describing  constraints. 

In  addition,  a  new  tuple  set  Valid  Beddown,  VB,  is  defined  in  order  to  describe  the 
added  summation  in  the  objective  function.  The  set  VB  =  {(/,  m)  |  H{i,  m)  =  1}  describes 
the  set  of  2-tuples  that  are  defined  only  if  Requirement  n  may  outload  at  POD  i  onto 
vehicle  Mode  m  on  any  day.  Detail  on  how  this  set  is  derived  is  given  in  (30)  below. 

II,  if  3  a  requirement  n  that  may  outload  at  POD  i  onto  a 

Mode  m  vehicle  on  any  day  (30) 

0<  otherwise 
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The  addition  of  the  MPBR  bears  similarity  to  the  GIBR  addition.  Both  attempt  to 
accomplish  the  same  effects  on  the  ITDM,  and  both  involve  a  single  addition  to  the 
objective  function  and  a  few  extra  constraints.  Again,  definition  of  the  beddown 
Mini/Max  decision  variables  Qim  need  not  be  defined  by  specific  POD  i.  They  can  also  be 
described  across  the  entire  theater  of  operations  including  all  PODs  in  a  TPFDD  utilized 
for  analysis.  In  other  words,  if  analysts  are  interested  in  minimizing  the  maximum  on 
beddowns  of  vehicle  Mode  m  by  Location  i,  the  variables  Qm  would  be  included  in  the 
model.  Further,  if  they  are  interested  in  doing  this  on  a  theater  wide  beddown  of  vehicle 
Mode  m,  the  variables  Qm  are  included.  This  describes  the  upper  bound  of  vehicles  of 
Mode  m  to  be  beddown  in  a  specific  theater  (all  PODs  in  TPFDD)  on  any  day.  Several 
more  variations  of  the  Mini/Max  Programming  decision  variables  exist,  such  as  Qm 
which  minimizes  the  upper  bound  on  the  number  of  vehicles  of  Mode  m.  Type  k  at  POD  i 
on  any  day,  and  Qmk,  which  minimizes  the  upper  bound  on  the  number  of  vehicles  of 
Mode  m.  Type  k  at  any  POD  on  any  day.  Note  that  as  the  decision  variables  are  indexed 
differently,  the  cost  parameter  d  must  be  indexed  differently  as  well.  For  the  purposes  of 
this  research  only  the  formulations  including  variables  Qim  and  Qm  and  cost  parameter  dm 
will  be  described.  The  Qim  defined  model  is  given  first.  Any  additions  made  to  the  ITDM 
is  outlined  and  described  in  Model  5  below.  See  Model  3  for  information  on  and 
formulation  of  the  ITDM. 
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Mini/Max  Programming  Beddown  Reduction  (MPBR) 


The  Mini/Max  Programming  Beddown  Reduetion  variation  of  the  ITDM  whieh 
describes  goal  variables  by  Location  i  and  Mode  m  is  formulated  below  in  Model  5. 
Tables  8  and  9  outline  new  parameters  and  decision  variables.  Additions  to  the  ITDM 
(31)-(33)  are  followed  by  an  asterisk. 

ITDM  With  MPBR  Addition  (by  Location  i,  Mode  m). 


Table  8.  MPBR  Parameters 


Parameter 

Description 

dm 

Cost  of  bedding  down  a  vehicle  of  Mode  m  at  any  POD 

Table  9.  MPBR  Decision  Variables 


Variable 

Description 

Qim 

Upper  bound  on  number  of  Mode  m  vehicles  to  be  beddown  at  POD  i  on 
any  day 

MPBR  Formulation. 

The  MPBR,  which  is  the  ITDM  model  with  Min/Max  additions,  is  shown  below 
in  Model  5. 
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(31) 


Minimize  2  +  X 

(i,J,m,k  ,v)  e  W  (n,i ,j,m,k  ,v)  e  LF 

y'  d  Q  * 

m^im 

{i,m)£VB 

Subject  to 

rd^  +q^ 


ZEZ  =  r,ij  V(n,  /,  j)  ^VR  (19) 

M  jj  v=ad„  +1 

^ijmkv  -^imv  \f{i,m,v)^VO  (20) 

%mk^ijmkv  -^jmv  \f {j ,m,V)  ^VU  (21) 

i  Km 

^  ynijmkv  -  ^ijmkv'^ijmkPk  V(/,  j,  M,  k,v)^W  (22) 

Nil, 

E  y(i,m,v)^VO,iyk^KJ  *  (32) 

ynijmkv  ^  0  V(n,  i,  j,  m,  k,  v)  &  VF  (23) 

a,„e{0}uZ"  \f(i,m)^VB  *  (33) 

^ijmkv  ^  {0}  ^  ^(b  2,  "J,  k,v)  (24) 


Model  5.  ITDM  With  MPBR  Addition  (By  Location  i,  Mode  m) 
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Addition  of  the  MPBR  presents  another  signifieant  improvement  on  the  beddown 
approximations  of  the  ITDM,  whieh  stem  direetly  from  its  vehiele  alloeation  solutions. 
First,  an  addition  to  the  objeetive  function  (31)  ensures  that  a  cost  is  assessed  whenever  a 
vehicle  of  Mode  m  is  beddown  at  a  POD.  When  this  cost  is  appropriately  set,  this  portion 
of  the  objective  function  ensures  a  smooth  mix  of  vehicle  allocations  across  a  given  time 
window.  The  objective  function  is  now  threefold.  It  attempts  to  simultaneously  minimize 
vehicle  usage  costs,  late  delivery  penalties,  and  maximum  beddown  costs  of  vehicle  Mode 
m.  We  will  see  later  on  how  manipulating  the  costs  bu  and  dm,  and  penalty  g  can  cause  the 
model  to  achieve  certain  desired  solutions. 

Two  sets  of  constraints  are  added  to  the  ITDM.  Constraints  at  (32)  ensure  that  the 
number  of  vehicles  required  to  make  all  deliveries  from  a  POD  i  is  less  than  the  maximum 
of  vehicles  at  a  specific  location.  When  these  constraints  are  combined  with  (30)  in  the 
objective  function,  the  model  attempts  to  minimize  the  maximum  number  of  vehicles  of 
Mode  m  at  POD  i  prescribed  in  a  solution,  resulting  in  a  reduced  beddown  approximation. 
Constraints  at  (33)  ensure  that  maximum  goal  variables  are  non-negative  and  integer,  and 
are  required. 

MPBR  Conclusion. 

The  MPBR  is  the  second  contribution  of  this  thesis.  When  combined  with  the 
initial  formulation  of  the  ITDM,  the  MPBR  provides  better  beddown  solutions  for  force 
flow  analysts  to  work  with.  Instead  of  solutions  accounting  only  for  daily  vehicle  costs, 
results  are  now  based  on  the  simultaneous  minimization  of  daily  costs  and  beddown  costs. 
The  model  now  attempts  to  give  a  minimal  vehicle  solution  based  on  two  cost  objectives 
as  well  as  the  reduction  of  late  deliveries.  Thus,  a  more  realistic  theater  distribution 
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modeling  affect  is  achieved.  While  the  MPBR  and  GIBR  have  different  formulations,  both 
achieve  the  same  effect.  For  this  reason,  only  MPBR  test  results  are  included  in  Chapter  4 
of  this  research. 

ITDM  With  MPBR  By  Location  and  Theater  (A  Multi-Ohjective  Approach) 

MPBR  Multi-Objective  Approach  Overview 

As  discussed  previously,  the  MPBR  decision  variables  Q  may  be  indexed  by 
specific  Location,  Qim,  or  by  the  entire  theater  of  interest,  Qm.  It  turns  out  that  both  of 
these  minimizations  can  be  achieved  simultaneously.  This  allows  force  flow  analysts  to 
set  limits  on  the  number  vehicles  of  a  specific  Mode  m  at  each  POD  i,  as  well  as  in 
theater,  (V/  e  I). 

As  with  previous  models,  all  original  aspects  of  the  ITDM  including  parameters, 
decision  variables,  set  defining  binary  functions,  basic  sets,  and  function  derived  tuple  sets 
remain  the  same.  In  fact,  several  of  these  aspects  are  utilized  when  defining  this 
formulation  of  the  MPBR  addition.  One  function  derived  tuple  set  utilized  in  the  MPBR  is 
the  set  of  Valid  Beddowns,  VB,  which  appears  in  the  objective  function.  Since  the 
variables  Qm  are  defined  across  all  PODs  i,  A  new  function  derived  tuple  set  is  created  for 
constraint  formulation.  The  set  VM  =  {(m,v)  |  L{m,v)  =  1}  describes  the  set  of  2-tuples 
that  are  defined  only  if  Requirement  n  may  outload  at  any  POD  in  theater  onto  vehicle 
Mode  m  on  day  v.  Detail  on  how  this  set  is  derived  is  given  in  (34)  below. 

fl,  if  3  a  requirement  n  that  may  outload  at  any  POD  in  theater  onto  a 

Mode  m  vehicle  on  Day  v  (34) 

O'  otherwise 

In  addition,  the  formulation  employs  two  of  the  ITDM  basic  sets.  The  basic  sets 
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utilized  in  the  MPBR  addition  are  the  set  of  all  vehieles  of  Mode  m  that  are  also  Type  k. 
Km,  and  the  set  of  all  vehieles  of  Mode  m,  M.  The  sets  M  and  VB  are  employed  in 
deseribing  summations  in  the  objective  function  while  the  sets  Km  and  Whelp  formulate 
necessary  constraints.  Two  costs  are  now  defined,  dim  and  dm,  which  are  indexed  by 
location  and  theater,  respectively.  This  Multi-Objective  approach  provides  an  example  of 
how  versatile  the  MPBR  can  be  in  terms  of  modeling  policy  driven  preferences.  Model 
formulation  including  both  the  MPBR  objective  by  Location  i  and  whole  theater  is  shown 
below  in  Model  6.  Tables  8  and  10  outline  the  parameters  and  decision  variables. 


Table  8.  MPBR  Parameters 


Parameter 

Description 

dm 

Cost  of  bedding  down  a  vehicle  of  Mode  m  Type  k  at  any  POD 

dim 

Cost  of  bedding  down  a  vehicle  of  Mode  m  Type  k  at  POD  i 

Table  10.  MPBR  Decision  Variables  (Multi-Objective 


Variable 

Description 

Qim 

Upper  bound  on  number  of  Mode  m  vehicles  to  be  beddown  at  POD  i  on 
any  day 

Qm 

Upper  bound  on  number  of  Mode  m  vehicles  to  be  beddown  in  theater  (all 
PODs)  on  any  day 

MPBR  Multi-Objective  Formulation. 

The  MPBR,  the  ITDM  model  with  Min/Max  additions  by  location  and  theater  is 
shown  below  in  Model  6. 
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Minimize  ^  + 

(i,J,m,k  ,v)  e  W  (n 

'y  d  Q  +'y  d  Q 

(i,m)eVB  M 

Subject  to 

rd^  +q^ 

Z  Siv-rdJy„,j„,^ 

•,i  ,J,m,k  ,v)  e  LF 

* 

n 

(35) 

3^ nijmkv  ^nij 

M-fj  v=ad„  -1-1 

V  (n,i,j)^VR 

(19) 

^ijmk  ^ijmkv  ^imv 

J 

y{i,m,v)  ^VO 

(20) 

^ijmk  ^ijmkv  ^  jmv 

i  Km 

y{J,m,v)^VU 

(21) 

3^ nijmkv  ~  ^ijmkv  ^ijmk  Pk 

(22) 

j:(i,j,m,k,v)GW 

yii,m,v)  ^VO  ,  (Vk  ^  KJ* 

(32) 

X  <o 

'^ijmkv  — 

i,j:(i,j,m,k,v)GW 

yim,v)  ^VM  ,  (yk  ^  K„j  * 

(36) 

y  nijmkv  —  ^ 

V(n,  i,  j,  m,  k,  v)  e  VF 

(23) 

a,„e{o}uz" 

V(/,  m)  e  VB  * 

(33) 

e„e{o}uz" 

V  m  eM  * 

(37) 

e{0}uZ^ 

k,v) 

(24) 

Model  6.  ITDM  With  MPBR  Addition  By  Location  and  Theater 

The  MPBR  formulation  with  multi-objective  addition  to  the  ITDM  presents  a 
technique  that  models  two  important  aspects  of  the  research  described  to  this  point.  With 


57 


the  addition  of  both  objectives  simultaneously  seen  in  (35),  and  the  new  constraints  at 
(32)-(33),  (36)-(37),  the  ITDM  provides  a  detailed  mixed  integer  programming  model  for 
force  flow  analysts  at  USTRANSCOM  to  utilize.  Objective  additions  at  (35)  ensure  vehicle 
solutions  are  minimized  by  mode  at  each  POD  and  by  mode  for  the  whole  theater.  The 
parameter  dim  is  indexed  by  Location  i  and  allows  the  model  to  define  differing  beddown 
costs  based  on  POD.  This  helps  when  a  user  desires  to  use  fewer  vehicles  of  a  certain 
Mode  m  at  a  specific  POD  i. 

Two  sets  of  constraints  are  added  to  the  ITDM.  Constraints  at  (32)  ensure  that  the 
number  of  vehicles  required  to  make  deliveries  in  a  TPFDD  by  location  is  less  than  the 
maximum  number  of  vehicles  at  that  location.  Constraints  at  (36)  ensure  that  the  number 
of  vehicles  required  to  make  all  deliveries  in  a  TPFDD  is  less  than  the  maximum  on 
vehicles  for  the  theater.  These  constraints  help  minimize  location  and  theater  wide 
beddown  approximations,  and  allow  analysts  to  set  limits  on  measures.  Constraints  at  (33), 
(37)  ensure  that  decision  variables  are  non-negative  and  integer. 

As  with  other  models  presented  in  this  research,  this  model  can  be  manipulated  to 
provide  certain  deliberate  solutions  based  on  current  policy  and  changing  objectives. 
Manipulation  can  be  achieved  by  adjusting  the  values  of  bu,  g,  and  dm  ■  The  ITDM  now 
attempts  to  reduce  beddown  by  location,  and  seeks  to  minimize  the  maximum  allocation  of 
vehicles  in  theater  as  described  in  the  MPBR  addition.  The  formulation  and  basic 
mathematics  of  this  cumulative  model  remain  unchanged  from  the  original  MPBR  model. 

It  simply  expands  upon  it.  The  goal  here  is  to  provide  a  tool  that  optimizes  several  realistic 
aspects  of  theater  distribution  simultaneously.  Different  variations  of  Min/Max  decision 
variables  could  be  easily  added  to  the  model,  depending  on  vehicle  needs  and  constraints. 
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Conclusion 


This  chapter  has  extensively  detailed  the  model  additions  developed  in  this 
researeh,  namely  the  GIBR,  the  MPBR,  and  the  Min/Max  Multi-Objeetive  approaeh.  The 
improvements  that  eaeh  addition  makes  on  the  ITDM  were  also  discussed.  Approximating 
measures  for  beddowns  are  earned  over  from  the  original  ITDM  and  remain  unehanged. 
These  measures  will  be  utilized  in  depth  in  Chapter  4  when  eomparing  original  solutions 
to  the  ITDM  with  improved  solutions.  The  next  ehapter  of  this  thesis  will  entail 
implementation  of  the  MPBR  model  addition  on  several  different  test  cases. 
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IV.  Implementation  and  Results 


Implementation. 

The  ITDM  developed  by  Hafich  was  implemented  using  both  Microsoft  Excel 
2007  and  the  optimization  software  LINGO  13(Lmdo  Systems  Inc,  2012).  This  was  done 
with  a  decision  Support  System  (DSS)  built  in  the  excel  domain.  The  DSS  was  organized 
such  that  a  user  uploads  a  TPFDD  and  enters  all  input  parameters  necessary  to  define  the 
model.  Once  a  TPFDD  is  selected  and  all  parameters  entered,  the  DSS  uses  Visual  Basic 
for  Applications  (VBA)  code  to  process  data  and  write  the  mixed  integer  programming 
model  in  the  LINGO  13  environment.  The  model  is  then  solved  by  LINGO  13  and 
solution  data  is  passed  back  to  the  Excel  environment  in  a  readable  format.  All  original 
elements  of  this  DSS  were  developed  by  Hafich  with  the  assistance  of  Dr.  Jeffery  Weir  of 
the  Air  Force  Institute  of  Technology.  This  version  of  the  DSS  is  relatively  unchanged 
when  solutions  to  the  ITDM  are  referenced  for  comparison  testing.  Readers  are 
encouraged  to  see  Appendix  I  for  ITDM  VBA  code  updates  utilized  in  this  research. 

To  test  whether  the  newly  developed  models  in  this  research  produce  better 
solutions,  VBA  code  additions  were  made  to  the  DSS  and  implemented.  These  updates 
align  directly  with  the  MPBR  mathematical  changes  made  to  the  ITDM.  Thus,  the 
process  of  obtaining  solutions  to  the  math  programming  model  via  Excel  VBA  2007  and 
LINGO  13  remained  the  same.  The  differences  were  the  formulation  of  the  model  in  the 
LINGO  interface,  and  the  resulting  solutions  from  this  formulation.  All  testing  was 
conducted  on  a  Lenovo  Think  Center  M58  computer  running  Windows  Vista  (Service 
Pack  2)  with  two  Intel  Celeron  2.6GHz  processors  and  4  GB  of  RAM. 
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Because  the  MPBR  increases  model  size  and  complexity,  it  often  increases  the 
time  to  find  an  optimal  solution.  This  is  because  the  Mini/Max  Programming  model  tends 
to  bounce  back  and  forth  between  solutions  within  the  branch  and  bound  process.  A 
relative  optimality  tolerance  was  set  to  encourage  faster  solutions,  and  the  solver  was  set 
to  search  for  solutions  within  5%  of  the  true  optimal  for  one  minute.  If  an  optimal 
solution  was  not  found  within  one  minute,  a  feasible  solution  within  5%  of  the  Linear 
Program  Relaxation  lower  bound  was  reported  as  globally  optimal.  For  consistency,  the 
same  relative  optimality  tolerance  was  used  when  testing  both  the  ITDM  and  MPBR. 
Other  settings  imposed  on  LINGO  13  for  this  chapter  are  available  for  review  in 
Appendix  A. 

Model  Testing. 

For  this  analysis,  ITDM  beddown  solutions  were  tested  and  compared  with 
MPBR  solutions  for  6  test  cases.  All  test  cases  were  notional.  The  first  four  test  cases 
involved  varying  vehicle  mode  and  type  constraints  within  the  DSS.  The  fifth  test  case 
was  carried  out  to  analyze  the  effect  of  widening  time  windows  in  a  notional  TPFDD, 
while  the  last  case  looked  at  equal  operating  costs  per  short  ton  and  possible  policy 
driven  solutions.  Most  solutions  of  the  ITDM  were  found  quickly,  while  the  majority  of 
MPBR  solutions  took  the  entire  minute  to  solve. 

For  each  test  case,  a  smaller  notional  TPFDD  was  used  and  solutions  compared. 
This  was  the  exact  TPFDD  and  data  used  as  an  example  in  the  internal  research  paper  by 
Longhorn  &  Kovich  (2012).  For  Test  Case  6  a  similar  but  larger  notional  TPFDD  was 
also  implemented  and  results  compared.  For  each  case,  information  regarding  beddown 
solutions  was  collected  and  reported.  This  beddown  information  was  taken  directly  from 
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the  nonzero  vehiele  alloeation  deeision  variables  Xijmkv-  The  smaller  TPFDD  used  is 
shown  in  Table  1 1  below. 


Table  11.  Smaller  Notional  TPFDD 


Requirement 

POD 

Destination 

Short  Tons 

EAD 

RDD 

1 

il 

.jl 

500 

2 

4 

2 

il 

.jl 

250 

3 

5 

3 

il 

jl 

750 

4 

6 

4 

il 

jl 

200 

5 

7 

5 

il 

jl 

100 

6 

8 

6 

il 

j2 

600 

2 

5 

7 

il 

j2 

400 

3 

6 

8 

il 

j2 

200 

4 

7 

9 

il 

j2 

300 

5 

8 

10 

il 

j2 

500 

6 

9 

11 

i2 

jl 

500 

4 

5 

12 

i2 

jl 

400 

5 

6 

13 

i2 

jl 

300 

6 

7 

14 

i2 

j2 

1000 

3 

5 

15 

i2 

j2 

200 

5 

7 

16 

i2 

j2 

500 

7 

9 

The  TPFDD  lists  16  movement  Requirements,  two  PODs,  and  two  Destinations. 
Note  that  the  Short  Tons  eolumn  gives  the  values,  the  BAD  eolumn  gives  the  ad„ 
values,  and  the  RDD  eolumn  gives  the  rdn  values.  Note  also  that  the  possible  delivery 
days,  ineluding  extension  days,  (i.e.  the  set  B)  ranges  between  Day  3  and  Day  10 
(Hafieh).  The  larger  TPFDD,  although  too  large  to  include  here,  is  very  similar  to  the 
smaller  one  but  contains  64  Requirements,  two  PODs,  and  two  Destinations. 

Some  constant  assumptions  were  made  throughout  for  consistency  in  model 
testing.  First,  it  was  assumed  that  Wijmk  =  1  in  all  cases.  That  is,  cycle  values  were  always 
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set  to  one,  meaning  eaeh  vehicle  could  make  a  single  pickup  and  delivery  per  day.  For 
example,  if  w_iiji_AiR_c-5  =  1,  then  it  is  possible  for  a  C-5  to  make  a  single  delivery  from 
il  to  jl,  return  to  il  the  same  day,  and  be  available  for  use  the  following  day.  In  addition, 
for  most  testing,  an  arbitrarily  large  bound  (1000)  was  set  on  outloading  and  unloading 
parameters  while  bounds  of  30  was  set  for  others.  Setting  these  bounds  directly  implies 
the  beddown  upper  bound  limits  at  each  Location  i.  All  of  this  information  was  clearly 
defined.  For  the  purpose  of  this  research,  problem  size  (the  number  of  variables  and 
constraints)  and  objective  function  values  were  considered  irrelevant  and  were  not 
recorded.  Additionally,  all  solutions  produced  were  100%  on  time  as  a  result  of  a  large 
late  penalty  value. 

In  all  tests,  it  was  assumed  that  requirements  arrived  at  the  POD  on  the  FAD 
stated  in  the  TPFDD  and  require  one  day  of  processing.  Thus,  for  each  requirement,  adn 
is  set  to  the  requirement’s  FAD  +  1  day.  Additionally,  every  requirement  was  given  a 
single  extension  day  within  all  test  cases.  That  is,  qdn  =  1  for  all  requirements  (Hafich). 

Deriving  Beddown  Solutions. 

Potential  beddown  solutions,  as  discussed  in  Chapter  II,  were  derived  from 
vehicle  allocation  variables  using  the  equation 


Beddown^^j,  —  max 

veV 


ijmkv 


V  J 


(25) 


All  results  reported  in  this  chapter  were  obtained  by  post-processing  solutions 
from  the  DSS  by  way  of  equation  (25). 
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Test  Case  1:  Comparison  Utilizing  All  Modes,  All  Types. 


As  with  all  test  cases,  the  first  test  case  utilized  the  same  TPFDD  from  the 
Longhorn  &  Kovich  paper  shown  in  Table  11.  Here,  all  Modes  (AIR,  RAIL,  ROAD),  and 
all  Types  (C-130,  C-17,  C-5,  HEMTT,  M1083,  M35,  DODX,  ITTX,  and  FTTX)  of 
vehicles  are  used  to  investigate  how  the  models  react.  The  penalty  per  day  per  late  short 
ton  was  set  to  g  =  10,000.  Daily  cost  and  payload  data  were  given  by  USTRANSCOM 
directly  and  are  notional.  The  payload,  cost,  and  unloading  parameters  used  are  shown  in 
Appendix  B.  After  the  ITDM  and  ITDM  with  MPBR  addition  were  tested  on  this  case, 
model  outputs  were  post-processed  for  beddown  approximations  and  these  solutions 
included  in  Tables  12  and  13  below. 


Table  12.  Beddowns  by  Vehicle  Type,  POD  For  Test  Case  l(Outload/UnloadlOOO) 


Model 

C-130 

c-17 

C-5 

HEMT 

M1083 

M35 

DODX 

ITTX 

ETTX 

TOTAL 

ITDM 

il 

0 

0 

0 

0 

0 

131 

0 

0 

0 

131 

i2 

0 

0 

0 

0 

0 

188 

0 

0 

0 

188 

MPBR 

il 

1 

1 

1 

1 

1 

1 

0 

1 

4 

11 

i2 

1 

2 

5 

1 

1 

0 

0 

1 

4 

15 

Table  13.  Beddowns  by  Vehicle  Type,  POD  For  Test  Case  1  (Outload/Unload  30) 


Model 

C-130 

C-17 

C-5 

HEMT 

M1083 

M35 

DODX 

ITTX 

ETTX 

TOTAL 

ITDM 

il 

30 

8 

0 

0 

0 

30 

0 

0 

0 

68 

i2 

30 

12 

0 

0 

0 

30 

0 

0 

0 

72 

MPBR 

il 

1 

1 

1 

1 

1 

1 

0 

1 

4 

11 

i2 

1 

2 

5 

1 

1 

0 

0 

1 

4 

15 

64 


Note  that  for  both  unload/outloud  bound  settings  (1000  and  30),  solutions  to  the 
ITDM  and  MPBR  are  drastieally  different  in  terms  of  total  vehicle  beddowns.  In  the  first 
case,  the  ITDM  produced  an  all  M35  solution  with  131  M35s  at  POD  il  and  188  M35s  at 
POD  i2.  Meanwhile,  the  MPBR  outputs  a  smooth  distribution  of  vehicles  with  a  total 
beddown  of  1 1  vehicles  at  POD  il  and  15  at  POD  i2,  presenting  a  significant 
improvement  over  the  ITDM.  Delivery  of  requirements  in  the  notional  TPFDD  is  spread 
more  efficiently  over  the  allowable  time  window.  For  this  test  case,  when  outload/unload 
parameters  are  constrained  at  30  as  shown  in  Table  13,  the  MPBR  still  produced 
significantly  reduced  beddown  results.  Notice  that  the  ITDM  reached  its  upper  bound  on 
both  C-130s  and  M35s.  These  unnecessarily  large  vehicle  beddowns  for  the  ITDM  result 
from  respective  vehicle  operating  costs  bu.  This  is  because  the  ITDM  selects  the  cheapest 
vehicles  in  terms  of  bk  and  attempts  to  move  as  many  requirements  as  possible  using 
these  types  of  vehicles.  Model  output  solutions  from  this  case  are  included  for  reference 
in  Figures  2  and  3  below. 
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63  M35(s) 

leaving  POD 

11 

for 

destination 

J1 

on 

day 

4  (ROAD) 

500.00  Short  Tons 

of  Movement 

1 

4.00  Short  Tons 

of  Movement 

2 

50  M35(s) 

leaving  POD 

11 

for 

destination 

J2 

on 

day 

4  (ROAD) 

400.00  Short  Tons 

of  Movement 

7 

31  M35(s) 

leaving  POD 

11 

for 

destination 

J1 

on 

day 

5  (ROAD) 

246.00  Short  Tons 

of  Movement 

2 

2.00  Short  Tons 

of  Movement 

3 

100  M35(s) 

leaving  POD 

11 

for 

destination 

J2 

on 

day 

5  (ROAD) 

600.00  Short  Tons 

of  Movement 

6 

200.00  Short  Tons 

of  Movement 

8 

94  M35(s) 

leaving  POD 

11 

for 

destination 

J1 

on 

day 

6  (ROAD) 

748.00  Short  Tons 

of  Movement 

3 

4.00  Short  Tons 

of  Movement 

4 

37  M35(s) 

leaving  POD 

11 

for 

destination 

J2 

on 

day 

6  (ROAD) 

296.00  Short  Tons 

of  Movement 

9 

37  M35(s) 

leaving  POD 

11 

for 

destination 

J1 

on 

day 

7  (ROAD) 

196.00  Short  Tons 

of  Movement 

4 

100.00  Short  Tons 

of  Movement 

5 

63  M35(s) 

leaving  POD 

11 

for 

destination 

J2 

on 

day 

7  (ROAD) 

4.00  Short  Tons 

of  Movement 

9 

500.00  Short  Tons 

of  Movement 

10 

63  M35(s) 

leaving  POD 

12 

for 

destination 

J1 

on 

day 

5  (ROAD) 

500.00  Short  Tons 

of  Movement 

11 

125  M35(s) 

leaving  POD 

12 

for 

destination 

J2 

on 

day 

5  (ROAD) 

1000.00  Short  Tons 

of  Movement 

14 

50  M35(s) 

leaving  POD 

12 

for 

destination 

J1 

on 

day 

6  (ROAD) 

400.00  Short  Tons 

of  Movement 

12 

25  M35(s) 

leaving  POD 

12 

for 

destination 

J2 

on 

day 

6  (ROAD) 

200.00  Short  Tons 

of  Movement 

15 

38  M35(s) 

leaving  POD 

12 

for 

destination 

J1 

on 

day 

7  (ROAD) 

300.00  Short  Tons 

of  Movement 

13 

63  M35(s) 

leaving  POD 

12 

for 

destination 

J2 

on 

day 

9  (ROAD) 

500.00  Short  Tons 

of  Movement 

16 

Figure  2.  ITDM  Case  1  Solution  Test  Case  1 
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1  C130(s) 

leaving  POD 

11 

for 

destination 

J2 

on 

day 

3  (AIR) 

12.00  Short  Tons 

of  Movement 

6 

1  C5(s) 

leaving  POD 

11 

for 

destination 

J1 

on 

day 

4  (AIR) 

50.00  Short  Tons 

of  Movement 

1 

10.00  Short  Tons 

of  Movement 

2 

1  C5(s) 

leaving  POD 

11 

for 

destination 

J2 

on 

day 

5  (AIR) 

59.00  Short  Tons 

of  Movement 

6 

1.00  Short  Tons 

of  Movement 

8 

1  C5(s) 

leaving  POD 

11 

for 

destination 

J1 

on 

day 

6  (AIR) 

60.00  Short  Tons 

of  Movement 

3 

1  C17(s) 

leaving  POD 

11 

for 

destination 

J2 

on 

day 

7  (AIR) 

4.00  Short  Tons 

of  Movement 

9 

31.00  Short  Tons 

of  Movement 

10 

1  C130(s) 

leaving  POD 

11 

for 

destination 

J2 

on 

day 

8  (AIR) 

12.00  Short  Tons 

of  Movement 

10 

2  C5(s) 

leaving  POD 

12 

for 

destination 

J2 

on 

day 

4  (AIR) 

120.00  Short  Tons 

of  Movement 

14 

2  C17(s) 

leaving  POD 

12 

for 

destination 

J1 

on 

day 

5  (AIR) 

50.00  Short  Tons 

of  Movement 

11 

1  C17(s) 

leaving  POD 

12 

for 

destination 

J2 

on 

day 

5  (AIR) 

35.00  Short  Tons 

of  Movement 

14 

1  C5(s) 

leaving  POD 

12 

for 

destination 

J2 

on 

day 

5  (AIR) 

60.00  Short  Tons 

of  Movement 

14 

1  C130(s) 

leaving  POD 

12 

for 

destination 

J1 

on 

day 

6  (AIR) 

12.00  Short  Tons 

of  Movement 

12 

1  C17(s) 

leaving  POD 

12 

for 

destination 

J1 

on 

day 

6  (AIR) 

35.00  Short  Tons 

of  Movement 

12 

1  C5(s) 

leaving  POD 

12 

for 

destination 

J1 

on 

day 

6  (AIR) 

60.00  Short  Tons 

of  Movement 

12 

1  C5(s) 

leaving  POD 

12 

for 

destination 

J2 

on 

day 

7  (AIR) 

60.00  Short  Tons 

of  Movement 

15 

5  C5(s) 

leaving  POD 

12 

for 

destination 

J2 

on 

day 

8  (AIR) 

300.00  Short  Tons 

of  Movement 

16 

1  C17(s) 

leaving  POD 

12 

for 

destination 

J2 

on 

day 

9  (AIR) 

35.00  Short  Tons 

of  Movement 

16 

3  FTTX(s) 

leaving  POD 

11 

for 

destination 

J2 

on 

day 

3  (RAIL) 

450.00  Short  Tons 

of  Movement 

6 

3  FTTX(s) 

leaving  POD 

11 

for 

destination 

J1 

on 

day 

4  (RAIL) 

450.00  Short  Tons 

of  Movement 

1 

1  FTTX{s) 

leaving  POD 

11 

for 

destination 

J2 

on 

day 

4  (RAIL) 

150.00  Short  Tons 

of  Movement 

7 

1  FTTX{s) 

leaving  POD 

11 

for 

destination 

J1 

on 

day 

5  (RAIL) 

60.00  Short  Tons 

of  Movement 

2 

90.00  Short  Tons 

of  Movement 

3 

1  ITTX(s) 

leaving  POD 

11 

for 

destination 

J1 

on 

day 

5  (RAIL) 

180.00  Short  Tons 

of  Movement 

2 

2  FTTX{s) 

leaving  POD 

11 

for 

destination 

J2 

on 

day 

5  (RAIL) 

63.00  Short  Tons 

of  Movement 

6 

237.00  Short  Tons 

of  Movement 

7 

4  FTTX{s) 

leaving  POD 

11 

for 

destination 

J1 

on 

day 

6  (RAIL) 

600.00  Short  Tons 

of  Movement 

3 

2  FTTX(s) 

leaving  POD 

11 

for 

destination 

J1 

on 

day 

7  (RAIL) 

200.00  Short  Tons 

of  Movement 

4 

100.00  Short  Tons 

of  Movement 

5 

Figure  2.  ITDM  Case  1  Solution  Test  Case  1  (cont.) 
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2  FTTX(s) 

leaving  POD 

11 

for 

destination 

J2 

on 

day 

7 

(RAIL) 

199.00  Short  Tons 

of  Movement 

8 

101.00  Short  Tons 

of  Movement 

9 

1  ITTX(s) 

leaving  POD 

11 

for 

destination 

J2 

on 

day 

8 

(RAIL) 

180.00  Short  Tons 

of  Movement 

9 

3  FTTX(s) 

leaving  POD 

11 

for 

destination 

J2 

on 

day 

9 

(RAIL) 

450.00  Short  Tons 

of  Movement 

10 

4  FTTX(s) 

leaving  POD 

12 

for 

destination 

J2 

on 

day 

4 

(RAIL) 

600.00  Short  Tons 

of  Movement 

14 

1  ITTX(s) 

leaving  POD 

12 

for 

destination 

J2 

on 

day 

4 

(RAIL) 

180.00  Short  Tons 

of  Movement 

14 

3  FTTX(s) 

leaving  POD 

12 

for 

destination 

J1 

on 

day 

5 

(RAIL) 

450.00  Short  Tons 

of  Movement 

11 

2  FTTX(s) 

leaving  POD 

12 

for 

destination 

J1 

on 

day 

6 

(RAIL) 

293.00  Short  Tons 

of  Movement 

12 

2  FTTX(s) 

leaving  POD 

12 

for 

destination 

J1 

on 

day 

7 

(RAIL) 

300.00  Short  Tons 

of  Movement 

13 

1  FTTX(s) 

leaving  POD 

12 

for 

destination 

J2 

on 

day 

7 

(RAIL) 

140.00  Short  Tons 

of  Movement 

15 

1  FTTX(s) 

leaving  POD 

12 

for 

destination 

J2 

on 

day 

8 

(RAIL) 

150.00  Short  Tons 

of  Movement 

16 

1  M35(s) 

leaving  POD 

11 

for 

destination 

J2 

on 

day 

3 

(ROAD) 

8.00  Short  Tons 

of  Movement 

6 

1  M1083(s) 

leaving  POD 

11 

for 

destination 

J2 

on 

day 

4 

(ROAD) 

5.00  Short  Tons 

of  Movement 

7 

1  M35(s) 

leaving  POD 

11 

for 

destination 

J2 

on 

day 

5 

(ROAD) 

8.00  Short  Tons 

of  Movement 

6 

1  M35(s) 

leaving  POD 

11 

for 

destination 

J2 

on 

day 

6 

(ROAD) 

8.00  Short  Tons 

of  Movement 

7 

1  M35(s) 

leaving  POD 

11 

for 

destination 

J2 

on 

day 

7 

(ROAD) 

8.00  Short  Tons 

of  Movement 

9 

1  HEMTT(s) 

leaving  POD 

11 

for 

destination 

J2 

on 

day 

8 

(ROAD) 

7.00  Short  Tons 

of  Movement 

9 

1  HEMTT(s) 

leaving  POD 

11 

for 

destination 

J2 

on 

day 

9 

(ROAD) 

7.00  Short  Tons 

of  Movement 

10 

1  M1083(s) 

leaving  POD 

12 

for 

destination 

J2 

on 

day 

5 

(ROAD) 

5.00  Short  Tons 

of  Movement 

14 

2  M1083(s) 

leaving  POD 

12 

for 

destination 

J2 

on 

day 

8 

(ROAD) 

10.00  Short  Tons 

of  Movement 

16 

1  HEMTT(s) 

leaving  POD 

12 

for 

destination 

J2 

on 

day 

9 

(ROAD) 

5.00  Short  Tons 

of  Movement 

16 

Figure  3.  ITDM  with  MPBR  Solution  Test  Case  1 


It  is  clear  from  the  MPBR  model  output  that  requirement  distribution  has  been 
smoothly  spread  over  the  time  window  for  eaeh  requirement,  ultimately  resulting  in  a 
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reduced  vehicle  solution.  Similar  smoothing  of  delivery  will  be  seen  throughout  this 
chapter  when  the  MPBR  addition  is  implemented  in  the  remaining  test  cases. 


Test  Case  2:  Comparison  Using  All  Modes,  Single  Type. 

The  second  test  case  involved  constraining  both  models  to  include  all  Modes 
(AIR,  RAIL,  ROAD)  but  only  a  single  Type  (C-I30,  HEMTT,  and  DODX)  within  each 
Mode.  The  penalty  per  day  per  late  short  ton  was  set  to  g  =  10,000.  The  same  daily  cost 
and  payload  data  used  for  Test  Case  I  is  used  here  and  will  remain  constant  until  Test 
Case  6  is  presented.  The  payload,  cost,  and  unloading  parameters  used  are  attached  in 
Appendix  C.  After  the  ITDM  and  ITDM  with  MPBR  addition  were  tested  on  this  case, 
model  outputs  were  post-processed  for  beddown  approximations  and  these  solutions 
included  in  Table  14  below. 


Table  14.  Beddowns  by  All  Mode,  Type,  POD  Test  Case  2  (Outload/Unload  1000) 


Model 

C-130 

HEMT 

DODX 

TOTAL 

ITDM 

il 

0 

156 

0 

156 

i2 

0 

214 

0 

214 

MPBR 

il 

0 

1 

4 

5 

i2 

1 

20 

3 

24 

When  the  ITDM  and  MPBR  are  tested  with  all  modes  by  a  single  vehicle  type 
with  arbitrarily  high  outloading  and  unloading  values  (1000),  beddown  approximations 
are  once  again  improved  greatly.  The  ITDM  proposes  an  all  Road  (HEMTT)  solution  as  a 
result  of  this  vehicle  being  the  least  expensive  option  in  terms  of  operating  costs  bk.  POD 
il  requires  156  total  vehicles  to  deliver  requirements  while  POD  i2  requires  a  total  of  214 
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vehicles.  When  compared  to  the  MPBR  beddown  approximation  of  5  total  vehicles  at 
POD  il  and  24  vehicles  at  POD  i2,  it  is  clear  that  this  beddown  solution  for  all  modes  and 
a  single  type  is  largely  reduced  from  that  of  the  ITDM.  Further,  delivery  of  requirements 
in  the  notional  TPFDD  is  spread  more  efficiently  over  the  allowable  time  window.  Model 
output  solutions  are  included  for  reference  in  Appendix  H. 

Test  Case  3:  Comparison  Using  a  Single  Mode. 

The  third  test  case  involved  constraining  both  models  to  include  a  single  Mode 
(AIR)  and  all  Types  (C-130,  C-17,  and  C-5)  within  Mode  AIR.  The  penalty  per  day  per 
late  short  ton  was  still  set  to  g  =  10,000.  Payload,  cost,  and  unloading  parameters  used  are 
shown  in  Appendix  D.  After  the  ITDM  and  ITDM  with  MPBR  addition  were  tested  on 
this  case,  model  outputs  were  post-processed  for  beddown  approximations  and  these 
solutions  included  in  Table  15  below. 


Table  15.  Beddowns  by  Single  Mode,  POD  For  Test  Case  3(Outload/Unload  1000) 


Model 

C-130 

C-17 

C-5 

TOTAL 

ITDM 

il 

100 

0 

0 

100 

i2 

126 

0 

0 

126 

MPBR 

il 

3 

3 

11 

17 

i2 

5 

8 

11 

24 

The  ITDM  proposes  an  all  C-130  solution  as  a  result  of  this  being  the  least 
expensive  option  in  terms  of  operating  costs  bk.  POD  il  requires  100  total  AIR  vehicles 
to  deliver  requirements  while  POD  i2  requires  a  total  of  126  AIR  vehicles.  When 
compared  to  the  MPBR  beddown  approximation  of  17  AIR  vehicles  at  POD  il  and  24 


AIR  vehicles  at  POD  i2,  it  can  be  inferred  that  this  Mode  AIR  beddown  solution  is  more 
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efficient  and  incurs  a  lower  overall  cost  than  that  of  the  ITDM.  To  illustrate  this  fact, 
Figures  4-5  and  Table  16  are  included  below.  Table  16  shows  total  vehicle  trips  (sum  of 
bars  in  Figures  4  and  5),  which  differs  from  the  beddown  approximations. 
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3  4  5  6  7  8  9 
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Figure  4.  ITDM  C-130  Allocations  By  Day  at  POD  il 


3  4  5  6  7  8  9 

Day 
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Figure  5.  MPBR  Mode  AIR  Allocations  By  Day  at  POD  il 
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Table  16.  Vehicle  Allocations  and  Cost  Information 


at  ] 

POD  il  For  Test  Case  3 

Model 

Type 

k 

Daily  Cost 

K 

Total  Vehicle 
Trips  V  V 

^ijmkv 

Total  Operation 
Cost 

^ijmkv 

ITDM 

C-130 

3 

317 

951 

C-17 

9 

0 

C-5 

16 

0 

MPBR 

C-130 

3 

5 

1,011 

C-17 

9 

4 

C-5 

16 

60 

It  remains  to  determine  whether  total  cost  in  terms  of  vehicle  trips  and  beddowns 
is  reduced  by  the  MPBR.  Comparison  of  Figures  4  and  5  reveals  that  spikes  in  vehicle 
allocations  is  greatly  reduced.  There  is  a  trade  off,  however,  in  terms  of  operating  costs 
evidenced  by  Table  28.  Although  total  operating  costs  increased  slightly  with  the  MPBR, 
its  solution  is  more  desirable  as  the  theater  wide  beddown  of  aircraft  was  cut  from  226  C- 
130s  with  ITDM  to  41  various  aircraft  as  seen  in  Table  15.  Note,  however,  that  the  cost  to 
beddown  a  C-130  is  almost  certainly  cheaper  than  this  cost  for  a  C-5  or  C-17.  Thus,  an  all 
C-130  beddown  solution  should  come  with  a  lower  cost  than  a  mixed  AIR  beddown 
approximation  equivalent  in  value.  As  with  Test  Cases  1  and  2,  delivery  of  requirements 
in  the  notional  TPFDD  is  spread  more  effectively  over  the  allowable  time  window. 
Similar  results  can  be  expected  for  testing  with  other  Modes  (ROAD,  RAIL).  Model 
outputs  are  included  for  reference  in  Appendix  H. 
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Test  Case  4:  Comparison  Using  a  Single  Mode,  Single  Type. 


The  fourth  test  ease  involved  eonstraining  both  models  to  inelude  a  single  Mode 
(AIR)  and  a  single  Type  (C-5)  within  Mode  AIR.  The  penalty  per  day  per  late  short  ton 
was  again  set  to  g  =  10,000.  The  payload,  eost,  and  unloading  parameters  used  are  shown 
in  Appendix  E.  After  the  ITDM  and  ITDM  with  MPBR  addition  were  tested  on  this  ease, 
model  outputs  were  post-proeessed  for  beddown  approximations  and  these  solutions  are 
ineluded  in  Table  17  below. 


Table  17.  Beddowns  by  Single  Type,  POD  For  Test  Case  4  (Outload/Unload  1000) 


Model 

C-5 

TOTAL 

THEATER 

ITDM 

il 

18 

18 

44 

i2 

26 

26 

MPBR 

il 

11 

11 

24 

i2 

13 

13 

These  results  present  perhaps  the  most  undeniable  validation  that  the  MPBR 
reduees  overall  theater  distribution  eost.  Solutions  of  previous  test  eases  showed  a 
mixture  of  vehieles  with  differing  daily  operating  eosts,  while  the  solution  to  Test  Case  4 
is  limited  to  Type  C-5  aireraft  having  eonstant  operating  eosts.  The  MPBR  produees  a  7 
vehiele  reduction  at  POD  il  while  it  shows  a  13  vehicle  reduction  at  POD  i2,  making  the 
theater  wide  beddown  improvement  20  C-5s  for  this  small  test  case.  Even  larger 
reductions  can  be  expected  for  larger  TPEDDs. 

To  illustrate  that  the  MPBR  allocates  vehicles  in  a  cheaper  and  smoother  manner. 
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figures  6  and  7  are  included  below.  The  “peanut  butter”  spread  effect  of  Mini/Max 
Programming  is  clearly  evident  with  investigation  of  these  figures. 


Figure  6.  C-5  Allocations  By  Day  at  POD  il 


1 

1 

1 II II  1 1. 

4  5  6  7  8  9 

Day 

mlTDM  mMGPBR 


Figure  7.  C-5  Allocations  By  Day  at  POD  i2 
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It  is  seen  from  the  figures  that  at  both  POD  il  and  i2,  the  MPBR  allocates  vehicles  more 
smoothly  over  the  time  window  than  the  ITDM  in  the  single  vehicle  type  case.  In 
addition,  total  C-5  trips  for  both  models  at  POD  il  is  equal  with  exactly  64,  while  total  C- 
5  trips  at  POD  i2  matches  for  the  two  models  as  well  with  5 1 .  This  indicates  that  total 
operating  costs  to  deliver  requirements  are  equivalent  at  Location  iias  well  as  i2. 
However,  cost  savings  come  from  the  significant  theater  wide  beddown  reduction  of  20 
C-5s  seen  in  Table  17.  Thus,  results  here  can  be  translated  directly  to  a  truly  lower  cost 
solution  in  terms  of  both  vehicle  operation  and  beddowns.  These  results  indicate  that  the 
MPBR  optimizes  cumulative  costs  when  a  single  type  is  utilized  within  the  math 
programming  models.  Similar  results  can  be  excepted  when  a  single  type  of  a  different 
mode  (ROAD,  RAIL)  is  tested.  Model  output  solutions  are  included  for  reference  in 
Appendix  H. 

Test  case  5:  Analysis  ofTPFDD  with  Wider  Time  Windows. 

It  still  remains  to  validate  that  when  the  MPBR  is  given  the  opportunity,  it  will 
spread  delivery  of  requirements  even  more  efficiently  over  a  wider  time  window.  That  is, 
the  wider  requirement  time  windows  are  within  a  TPFDD,  the  fewer  the  number  of 
vehicles  required  to  deliver  its  requirements  will  be.  This  is  a  feat  that  the  ITDM  is  not 
able  to  accomplish.  Up  to  this  point,  testing  has  been  done  on  a  TPFDD  with  relatively 
narrow  time  windows  averaging  3  days.  Test  Case  5  involves  comparing  solutions  of  the 
MPBR  using  the  original  notional  TPFDD  to  solutions  using  the  same  TPFDD  with 
wider  requirement  delivery  time  windows.  The  only  changes  made  to  the  updated 
TPFDD  are  increases  of  7  in  {adn+1  to  rJ„+  qdn  )  for  each  respective  Time  Window. 
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This  is  done  by  increasing  the  value  of  the  Required  Delivery  Date  rdn  by  7  throughout 
the  TPFDD.  This  TPFDD  is  included  in  Table  18  below. 


Table  18.  Smaller  Notional  TPFDD 
With  Large  Time  Windows 


Requirement 

POD 

Destination 

Short  Tons 

EAD 

RDD 

1 

il 

.jl 

500 

2 

11 

2 

il 

jl 

250 

3 

12 

3 

il 

jl 

750 

4 

13 

4 

il 

jl 

200 

5 

14 

5 

il 

jl 

100 

6 

15 

6 

il 

j2 

600 

2 

12 

7 

il 

j2 

400 

3 

13 

8 

il 

j2 

200 

4 

14 

9 

il 

j2 

300 

5 

15 

10 

il 

j2 

500 

6 

16 

11 

i2 

jl 

500 

4 

12 

12 

i2 

jl 

400 

5 

13 

13 

i2 

jl 

300 

6 

14 

14 

i2 

j2 

1000 

3 

12 

15 

i2 

j2 

200 

5 

14 

16 

i2 

j2 

500 

7 

16 

As  with  all  test  cases,  the  fifth  test  case  utilized  the  same  TPFDD  from  the 
Longhorn  &  Kovich  paper  shown  in  Table  1 1.  It  also  implements  the  modified  version  of 
the  notional  TPFDD  at  Table  15.  Here,  all  Modes  (AIR,  RAIL,  ROAD),  and  all  Types 
(C-130,  C-17,  C-5,  HEMTT,  M1083,  M35,  DODX,  ITTX,  and  FTTX)  of  vehicles  are 
used  to  investigate  how  the  MPBR  model  reacts.  The  penalty  per  day  per  late  short  ton 
was  set  to  g  =  10,000  once  again.  Payload,  cost,  and  unloading  parameters  used  are 


76 


shown  in  Appendix  F.  After  the  ITDM  and  ITDM  with  MPBR  addition  were  tested  on 
this  ease,  model  outputs  were  post-proeessed  for  beddown  approximations  and  these 
solutions  ineluded  in  Table  19  below. 


Table  19.  Beddowns  by  Vehicle  Type,  POD-Test  Case  5  (Outload/Unload  1000) 


Model 

C-130 

C-17 

C-5 

HEMT 

M1083 

M35 

DODX 

ITTX 

ETTX 

TOTAL 

MPBR 

Narrow 

Window 

il 

1 

1 

1 

1 

1 

1 

0 

1 

4 

11 

i2 

1 

2 

5 

1 

1 

0 

0 

1 

4 

15 

MPBR 

Wide 

Window 

il 

0 

0 

1 

0 

0 

0 

0 

0 

2 

3 

i2 

0 

0 

0 

0 

0 

1 

0 

0 

2 

3 

Note  that  the  Narrow  Window  Solution  matehes  the  MPBR  solution  to  Test  Case 
1 ,  as  nothing  changed  for  the  input  parameters  and  within  the  TPFDD  used.  However,  the 
Wide  Window  Solution  requires  significantly  fewer  vehicles.  All  requirements  within  the 
TPFDD  can  be  moved  with  6  vehicles  in  theater,  while  the  narrow  window  solution 
requires  26.  This  follows  intuitively,  as  fewer  vehicles  should  be  required  when  more 
time  is  allotted  to  distribute  requirements.  Thus,  the  MPBR  model  formulation  acts  as 
desired,  and  tends  to  spread  delivery  as  much  as  possible  over  a  time  window.  This 
presents  a  highly  desirable  result  for  force  flow  analysis  at  USTRANSCOM  because  in 
most  realistic,  operational  TPFDDS,  time  windows  are  much  wider  than  those  contained 
in  the  notional  TPFDD.  Model  output  solutions  are  included  for  reference  in  Appendix  H. 
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Test  Case  6:  Investigating  Equal  Operating  Costs  Per  Short  Ton 

Test  Cases  1-5  involved  setting  daily  operating  eosts  bk  to  USTRANSCOM 
provided  preferenees.  In  every  ease,  when  the  ITDM  was  tested,  there  was  a  elear  lowest 
eost  vehicle  for  the  ITDM  mixed  integer  program  to  select.  It  is  interesting  to  point  out 
that  the  ITDM  selects  vehicles  based  on  lowest  operating  cost  per  short  ton,  not  just 
lowest  operating  cost.  That  is,  the  ratio  of  operating  cost  to  average  pay  load,  bi/pk,  for  a 
Type  k  vehicle  actually  determines  whether  a  vehicle  is  the  cheapest  within  the  model. 
Since  all  ITDM  test  cases  prescribed  large  numbers  of  vehicles  with  characteristic  lowest 
cost  per  short  ton,  it  remains  to  investigate  the  effect  of  setting  bk/pk  constant  for  all 
vehicles  k. 

The  final  test  case  looked  at  the  effects  of  standardizing  operating  costs  per  short 
ton  of  payload.  The  results  of  the  ITDM  were  compiled  and  compared  with  MPBR 
findings.  First,  all  Modes  (AIR,  RAIL,  ROAD),  and  all  Types  (C-130,  C-17,  C-5, 
HEMTT,  M1083,  M35,  DODX,  ITTX,  and  FTTX)  of  vehicles  were  used  to  investigate 
how  both  models  react  with  the  larger  notional  TPFDD.  Both  models  are  then  constrained 
to  include  a  single  Mode  (AIR)  and  all  Types  (C-130,  C-17,  and  C-5)  within  Mode  AIR, 
and  tested  on  the  smaller  notional  TPFDD.  The  penalty  per  day  per  late  short  ton  was 
again  set  to  g  =  10,000.  The  payload  and  cost  parameters  are  included  in  Table  20  below. 
Unloading  and  outloading  parameters  are  shown  in  Appendix  F  for  reference.  After  the 
ITDM  and  ITDM  with  MPBR  addition  were  tested  on  both  the  smaller  and  larger 
TPFDD,  model  outputs  were  post-processed  for  beddown  approximations  and  these 
solutions  included  in  Tables  21  and  22  below. 
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Table  20.  Vehicle  Parameters  for  Test  Case  6 


Type 

Average  Payload 

Pk 

Daily  Cost 

K 

Cost  Per  Ton 

bk/pk 

Beddown  Cost 

dm 

C-130 

12 

8.4 

0.7 

5 

C-17 

35 

24.5 

0.7 

5 

C-5 

3^60 

42 

0.7 

5 

HEMTT 

35  7 

rn> 

4.9 

0.7 

5 

M1083 

3.5 

0.7 

5 

M35 

UU 

5.6 

0.7 

5 

DODX 

o 

o 

140 

0.7 

5 

FTTX 

105 

0.7 

5 

ITTX 

uO 

5180 

_ _ 

126 

0.7 

5 

f 


8 


Table  21.  Beddo\|ii  by  Vehicle  Type,  POD  For  Test  Case  1  (Larger  TPFDD) 


Model 

C-130 

f! 

1-17 

C-5 

HEMT 

M1083 

M35 

DODX 

ITTX 

ETTX 

TOTAL 

ITDM 

11 

DC 

30 

) 

1 

il 

0 

0 

1 

1 

1 

13 

10 

3 

29 

i2 

0 

y 

i 

2 

1 

1 

0 

15 

3 

8 

32 

MPBR 

il 

0 

0 

0 

0 

0 

12 

5 

4 

21 

i2 

1 

1 

0 

0 

0 

15 

6 

4 

31 

180 


Table  22.  Beddowns  by  Vehicle  Type,  POD  For  Test  Case  1  (Smaller  TPFDD) 


Model 

C-130 

C-17 

C-5 

TOTAL 

ITDM 

il 

25 

15 

17 

57 

i2 

25 

20 

5 

50 

MPBR 

il 

3 

3 

11 

17 

i2 

5 

8 

11 

24 
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Table  21  shows  results  from  testing  a  larger  TPFDD  eontaining  64  requirements. 
Setting  the  ratio  of  operating  cost  per  short  ton  to  a  constant  value  across  all  vehicle 
Types  k  clearly  improves  ITDM  beddown  approximations  in  terms  of  total  vehicles. 

Note,  however,  that  this  solution  prescribes  primarily  RAIL  vehicles,  which  realistically 
can  carry  a  high  initial  beddown  cost.  In  fact,  the  MPBR  prescribes  a  largely  RAIL 
solution  as  well.  It  seems  here  that  when  bk/pk  is  constant  for  all  Types  k  both  the  ITDM 
and  MPBR  tend  to  select  vehicles  capable  of  delivering  higher  payloads.  The  MPBR  still 
requires  fewer  vehicles  to  deliver  requirements  with  a  total  in  theater  beddown  of  52 
vehicles.  This  is  a  9  vehicle  improvement  over  the  ITDM  solution  of  61  vehicles.  Thus, 
even  with  constant  costs  per  short  ton  the  Mini/Max  approximation  provides  a  reduced 
beddown  approximation. 

More  significant  improvements  are  seen  in  the  single  mode  (AIR)  test.  Testing 
was  performed  on  the  smaller  notional  TPFDD  to  show  that  the  MPBR  present  a  reduced 
mix  of  AIR  vehicles  across  the  board.  Although  the  ITDM  does  select  a  mixture  of  the 
vehicles  available  (50  C-130s,  35  C-17s,  22  C-5s)  instead  of  strictly  C-130s,  the  numbers 
are  unnecessarily  high.  This  fact  is  validated  by  the  MPBR  solution  of  AIR  vehicles  (8  C- 
130s,  1 1  C-17s,  22  C-5s).  Since  model  input  remained  unchanged,  this  solution  is  exactly 
the  same  as  the  MPBR  results  from  Test  Case  3.  Thus,  although  setting  operating  costs 
per  short  ton  equal  does  improve  upon  single  mode  beddown  approximations  given  by 
the  ITDM,  it  does  not  improve  single  mode  solutions  from  the  MPBR.  Additionally, 
because  of  vehicle  beddown  reductions,  the  MPBR  provides  a  better  solution  than  the 
ITDM  here. 
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Verification  and  Validation 


Any  time  one  builds  a  new  model,  it  must  be  verified  and  validated  before 
analysis  ean  be  eonsidered  eomplete.  Verifieation  ensures  that  the  model  is  being  built 
the  right  way.  Validation  ensures  that  the  right  model  is  being  built.  These  teehniques  are 
diseussed  for  the  two  models  presented  and  tested  in  Chapter  4  of  this  researeh. 

Verification. 

The  ITDM  is  easily  verified  by  the  results  from  Test  Case  1.  First,  Sinee  the 
ITDM  formulation  and  aeeompanying  DSS  were  taken  direetly  from  researeh  presented 
by  Hafieh  (2013),  its  solutions  within  this  researeh  are  exaetly  the  same  as  solutions 
given  by  Hafieh,  assuming  identieal  input  parameters.  Furthermore,  before  testing  was 
done,  one  of  the  ITDM  solutions  given  by  Hafieh  was  identieally  reprodueed  with  the 
model  used  for  this  researeh.  Thus,  the  ITDM  is  verified.  The  MPBR  is  verified  by  the 
faet  that  in  Test  Case  4,  the  exaet  same  number  of  total  vehiele  alloeations  was  output.  It 
is  also  verified  in  seeing  that  vehiele  alloeations  are  similar  to  the  ITDM  in  Test  Case  6, 
although  signifieantly  redueed  due  to  the  Min/Max  formulation.  Additionally,  the  ITDM 
maintains  its  mixed  integer  programming  formulation  with  inelusion  of  the  MPBR.  It 
also  produeed  feasible  solutions  to  this  model  in  LINGO  13  for  all  input  settings 
investigated  in  Test  Cases  1-6. 

Validation. 

Test  Case  1  offers  one  reason  why  the  ITDM  eannot  be  the  best  model  for 
approximating  beddowns  in  theater.  An  unneeessarily  large  number  of  vehieles  is 
preseribed  to  deliver  the  16  requirements.  With  the  MPBR  addition,  the  model  avoids 
alloeation  of  unneeessary  vehieles,  resulting  in  a  greatly  redueed  beddown  solution.  Test 
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Case  1  validates  MPBR  effectiveness  when  all  modes  and  types  are  utilized,  while  Test 
Cases  2-4  validate  other  mode/type  combinations.  The  ITDM  fails  to  provide  a  “peanut 
butter  spread”  of  vehicle  allocations  throughout  a  given  time  window.  The  MPBR 
corrects  this  issue.  The  most  notable  validation  of  this  fact  is  seen  in  Test  Case  5,  where 
fewer  vehicles  were  allocated  to  deliver  requirements  with  wider  time  windows. 

The  fact  that  ITDM  solutions  are  affected  so  much  by  variations  in  operating 
costs  demonstrates  that  it  is  not  entirely  useful  as  formulated.  As  operating  costs  are 
greater  for  larger  vehicles,  and  vary  with  USTRANSCOM  policy,  an  effective  model  is 
one  that  is  not  overly  sensitive  to  changes  in  these  costs.  The  MPBR  addresses  this,  as  its 
solutions  react  to  beddown  costs  as  well  as  operating  costs.  This  fact  is  validated  by 
results  from  Test  case  6. 
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V.  Conclusions  and  Future  Research 


Conclusions 

The  three  model  variations  presented  in  this  researeh  involved  additions  to  the 
objective  function  and  constraints  of  the  ITDM.  First,  the  GIBR  added  a  penalty  for 
increasing  beddown  size  throughout  a  time  window.  It  also  constrained  the  number  of 
vehicles  used  on  a  certain  day  to  being  no  more  than  the  number  of  vehicles  used  on  the 
previous  day.  Next,  the  MPBR  applied  a  penalty  to  the  beddown  of  each  mode  of  vehicles. 
It  also  constrained  the  number  of  vehicle  allocations  to  being  less  than  a  certain  number, 
which  was  minimized  in  the  objective  function.  In  terms  of  beddown  approximating,  the 
GIBR  and  MPBR  produce  the  same  results.  However,  it  is  unclear  whether  they  will 
produce  the  exact  same  vehicle  allocation  solutions.  Finally,  the  Multi-Objective  Approach 
provided  is  just  one  of  many  possible  existing  variations  to  the  MPBR.  It  focuses  on 
minimizing  beddowns  at  each  POD  as  well  as  the  total  in  theater  beddown. 

The  GIBR  and  the  MPBR  improve  upon  the  ITDM  making  it  a  much  more  useful 
tool  for  theater  distribution  analysis.  As  the  ITDM,  the  baseline  model,  was  tested  and 
analyzed,  it  became  clear  that  beddown  approximations  would  be  much  too  large  to  be 
realistically  feasible.  To  address  this,  the  GIBR  was  developed,  which  although  never 
tested,  provided  an  improved  formulation  of  the  ITDM.  This  formulation  was  expected  to 
spread  requirement  distribution  more  efficiently  over  a  given  time  window.  However,  the 
GIBR  created  excess  decision  variables  and  constraints,  as  its  accompanying  decision 
variables  were  indexed  by  day.  Thus,  the  MPBR  was  developed,  which  provided  the  same 
modeling  effects  of  the  GIBR,  but  decision  variables  that  were  not  defined  by  day.  This 
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reduces  the  problem  size  and  increases  solution  speed  over  the  GIBR. 


The  MPBR  gives  force  flow  analysts  a  much  better  approximation  of  vehicles 
needed  for  theater  distribution.  In  terms  of  problem  size,  the  MPBR  increases  the  size  of  the 
ITDM,  but  not  significantly.  Solutions  take  longer  to  produce,  but  a  relative  optimality 
tolerance  can  easily  be  set  to  encourage  faster  results.  Additionally,  the  integrality  property 
of  decision  variables  can  be  relaxed  to  speed  up  the  model.  The  ITDM  forces  unpredictable 
spikes  in  vehicle  allocations,  even  when  equivalent  operating  costs  per  short  ton  are  used  as 
input.  When  applied  to  the  ITDM,  the  MPBR  uses  an  iterative  process  to  take  these  vehicle 
allocation  spikes  on  certain  days  and  transfer  deliveries  to  other  days  where  fewer 
allocations  were  made.  This  is  because  a  onetime  penalty  is  applied  to  the  maximum 
number  of  vehicles  of  each  mode  prescribed  to  distribute  requirements.  Thus,  the  model 
attempts  to  minimize  high  vehicle  allocations.  The  result  is  a  much  smoother  and  more 
desirable  delivery  schedule,  as  beddown  approximations  are  greatly  reduced. 

The  ITDM  with  MPBR  is  able  to  find  feasible  vehicle  mixtures  that  minimize 
operational  cost,  late  deliveries,  and  beddown  size  simultaneously.  Since  costs  are  user 
defined,  solutions  can  be  vectored  toward  a  vehicle  mixture  that  aligns  with  current  policy. 
Furthermore,  because  various  cost  settings  drive  different  vehicle  mixtures,  results  may  be 
sensitive  to  alternate  optimal  solutions.  Thus,  one  beddown  solution  could  be  desired  over 
another  even  though  both  are  reported  as  optimal  cost  wise. 

The  MPBR  has  the  potential  to  provide  reduced  vehicle  mixtures  when  post¬ 
processing  results  of  the  ITDM  and  analyzing  possible  beddowns,  which  can  result  in  lower 
cost  starting  points  in  terms  of  vehicles  to  support  theater  distribution.  Through  this 

addition  to  the  ITDM  and  associated  Decision  Support  System,  force  flow  analysts  at 
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USTRANSCOM  possess  a  highly  flexible  tool  to  assist  in  theater  distribution  analysis.  The 
MPBR  portion  of  the  formulation  ean  be  easily  altered  based  on  preferenees  and  existing 
policy.  This  use  of  the  updated  ITDM  to  model  theater  distribution  and  estimate  beddowns 
has  the  potential  to  save  increasingly  valuable  vehicle  resources  and  DOD  funds. 

Future  Research 

There  are  many  areas  in  which  this  research  can  be  furthered  to  improve  upon 
results,  organize  testing  procedures,  and  model  operational  realities.  The  greatest  potential 
for  future  research  is  to  investigate  how  variations  in  the  weighting  of  costs  within  the 
objective  function  affect  solutions.  Through  the  research  process,  it  became  clear  that 
model  output  was  highly  sensitive  to  user  defined  costs  including  daily  operating  cost,  the 
late  delivery  penalty,  and  beddown  cost.  Although  solutions  were  somewhat  sensitive  to 
the  operational  costs  and  the  late  penalty,  they  were  most  affected  by  variations  in  the 
beddown  cost.  This  was  evident  when  multiple  mode  and  type  mixtures  were  used  as 
input.  For  example,  in  running  two  tests  with  the  exact  same  model  inputs,  even  a  small 
change  in  the  beddown  penalty  could  produce  significantly  different  vehicle  allocations. 
Research  into  the  effects  of  varying  the  different  model  costs,  and  organizing  these  effects 
in  some  standardized  format  should  prove  fruitful.  It  would  also  allow  force  flow  analysts 
to  predict  model  output  based  on  model  input,  which  could  save  significant  time  when  a 
specific  vehicle  solution  is  needed. 

An  investigation  of  how  vehicle  batching  affects  solutions  of  the  MPBR 
formulation  should  prove  to  be  practical  research.  Operational  realities  and  current  policy 
often  requires  vehicles  to  be  deployed  into  theater  in  batches.  This  is  because  vehicles 
typically  relocate  or  deploy  as  a  unit  containing  a  set  number  of  this  type  of  vehicle.  For 
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example,  trucks  usually  move  into  theater  with  the  United  States  Army  (USA)  as  a 
Battalion  or  Company,  while  aircraft  typically  move  with  the  United  States  Air  Force 
(USAF)  as  a  wing  or  squadron.  This  would  constrain  the  beddown  of  vehicles  of  a  certain 
type  to  be  a  multiple  of  the  number  of  vehicles  contained  in  that  unit  type’s  typical  batch 
size.  Investigating  this  issue  by  coordinating  with  USTRANSCOM  to  determine  batch  size 
should  add  to  the  realistic  effects  of  the  research. 

In  addition,  variations  in  the  formulation  of  the  MPBR  such  as  the  Multi-Objective 
Approach  should  be  tested.  It  should  be  determined  if  several  Mini/Max  beddown 
objectives  can  be  achieved  at  the  same  time,  as  higher  leadership  may  have  multiple 
preferences  for  vehicle  beddown  numbers. 

An  exploration  of  how  the  General  Integer  Beddown  Reduction  formulation 
presented  in  this  thesis  can  affect  solutions  would  provide  insight  into  whether  or  not  this 
beddown  improvement  technique  is  better  than  the  MPBR.  Although  the  same  modeling 
effect  is  expected  in  terms  of  beddown  approximations,  the  GIBR  may  provide  vehicle 
allocation  solutions  throughout  a  given  time  window  that  are  better  in  terms  of  operational 
costs.  This  would  involve  updating  the  VBA  code  within  the  ITDM  DSS  to  match  the 
modeling  changes  of  the  GIBR. 

Lastly,  as  mentioned  in  the  research  by  Hafich,  further  research  into  defining  cycle 
values  should  be  conducted.  Instead  of  relying  on  input  from  a  user  to  determine  feasible 
cycles,  a  tool  could  be  developed  that  accounts  for  operational  capabilities  such  as  vehicle 
speeds  and  outload/unload  times,  and  returns  a  specific  cycle  value.  This  pre-processing 
result  could  then  be  used  as  input  before  implementing  the  DSS. 
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Appendix  A.  LINGO  13  Settings  File  Contents 


The  LINGO. CNF  file  eontains  settings  whieh  have  been  changed  from  their 
default  values  within  LINGO  13.  The  contents  of  the  LINGO. CNF  file  as  utilized  in 
this  thesis  appear  below  (Hafich). 


Lingo  CNF  info: 

!  LINGO  Custom  Configuration  Data: 
MXMEMB=  25000 
ABSINT=  O.lOOOOOOOE-11 
IPTOLR=  0.50000000E-01 
TIM2RL=  120 
LINLEN=  150 
DUALCO=  0 
PRECIS=  12 
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Appendix  B.  Additional  Model  Inputs  for  Test  Case  1 


Table  23.  Vehicle  Parameters  for  Test  Case  1 


Type 

Average  Payload 

Pk 

Daily  Cost 

K 

Beddown  Cost 

dm 

C-130 

6 

3 

10 

C-17 

12 

9 

10 

C-5 

35 

16 

10 

HEMTT 

35  7 

1 

10 

M1083 

1 

10 

M35 

uu 

1 

10 

DODX 

^00 

60 

10 

FTTX 

42 

10 

ITTX 

5180 

_ cm _ 

52 

10 

Table  24.  Qodoading  Parameters  for  Test  Case  1 


1! 

IT 

Mode 

Outload  Capacity 

^imv 

li 

i 

Air 

1000 

1* 

m 

Road 

1000 

H 

Rail 

1000 

Air 

1000 

hP-i 

uoo 

OD 

Road 

1000 

Rail 

1000 

*Note,  for  each  POD/Mode^fftQr,  t 


j _ I _ I 

he  outload  capacity  is  assumed  constant  for  all  v 


Table  25.  Unloading  Parameters  for  Test  Case  1 


Destination 

Mode 

Unload  Capacity 

^  jmv 

Air 

1000 

Road 

1000 

Rail 

1000 

J2 

Air 

1000 

Road 

1000 

Rail 

1000 

*Note,  for  each  Destination/Mode  pair,  the  unload  capacity  is  constant  for  all  v. 
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Appendix  C.  Additional  Model  Inputs  for  Test  Case  2 


Table  26.  Vehicle  Parameters  for  Test  Case  2 


Type 

Average  Payload 

Pk 

Daily  Cost 

K 

Beddown  Cost 

dm 

C-130 

12 

3 

10 

HEMTT 

7 

1 

10 

DODX 

^^00 

60 

10 

§§ 

Table  27.  Outloading  Parameters  for  Test  Cases  2 


*Note,  for  each  POD/M 


P0D 


il 

-TgD- 


m- 


Mode 


Air 


Road 


Rail 


Air 


Road 


Rail 


od^gj^ir,  the  out 


Outload  Capacity 


1000 


1000 


1000 


1000 


1000 


1000 


oad  capacity  is  assumed  constant  for  all 


Days  V. 


Table  28.  Unloading  Parameters  for  Test  Case  2 


Destination 

Mode 

Unload  Capacity 

^  jmv 

Air 

1000 

Road 

1000 

Rail 

1000 

Air 

1000 

Road 

1000 

Rail 

1000 

*Note,  for  each  Destination/Mode  pair,  the  unload  capacity  is  assumed  constant  for 
all  Days  v. 


89 


Appendix  D.  Additional  Model  Inputs  for  Test  Case  3 


Table  29.  Vehicle  Parameters  for  Test  Case  3 


Type 

Average  Payload 

Pk 

Daily  Cost 

K 

Beddown  Cost 

dm 

C130 

12 

3 

5 

C-17 

35 

9 

5 

C-5 

3^60 

16 

5 

35 

60 

Table  30.  Oifiloading  Parameters  for  Test  Cases  3 


P@D 

Mode 

Outload  Capacity 

^imv 

Air 

1000 

7^ 

Air 

1000 

*Note,  for  each  POD/Mode  p^ggthe  outload  capacity  is  assumed  constant  for  all  days  v 

150 

1§§ 

180 

180 

Table  31.  Unloading  Parameters  for  Test  Case  3 


Destination 

Mode 

Unload  Capacity 

^  jmv 

Jl 

Air 

1000 

J2 

Air 

1000 

*Note,  for  each  Destination/Mode  pair,  the  unload  capacity  is  assumed  constant  for 
all  days  v . 
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Appendix  E.  Additional  Model  Inputs  for  Test  Case  4 


Table  32.  Vehicle  Parameters  for  Test  Case  4 


Type 

Average  Payload 

Pk 

Daily  Cost 

K 

Beddown  Cost 

dm 

C-5 

60 

16 

5 

635 


60 


Table  33.  Outloading  Parameters  for  Test  Cases  4 


POD 

5 

Mode 

Outload  Capacity 

^imv 

il 

Air 

1000 

a 

Air 

1000 

*Note,  for  each  POD/M 


od|^ir,  the  out: 


oad  capacity  is  assumed  constant  for  all 


days  V. 


150 


180 

Table  34.  Unloading  Parameters  for  Test  Case  4 


Destination 

Mode 

Unload  Capacity 

^  jmv 

Air 

1000 

J2 

Air 

1000 

*Note,  for  each  Destination/Mode  pair,  the  unload  capacity  is  assumed  constant  for 
all  days  v. 
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Appendix  F.  Additional  Model  Inputs  for  Test  Case  5 


Table  35.  Vehicle  Parameters  for  Test  Case  5 
TPFDD  With  Wide  Time  Windows 


Type 

Average  Payload 

Pk 

Daily  Cost 

K 

Beddown  Cost 

dm 

C130 

12 

3 

25 

C-17 

35 

9 

25 

C-5 

5^60 

16 

25 

HEMTT 

35  7 

1 

25 

M1083 

1 

25 

M35 

uu 

nh  ^ 

1 

25 

DODX 

^00 

60 

25 

FTTX 

ZO^DEO 

42 

25 

ITTX 

uO 

5180 

_ cm _ 

52 

25 

8 

200 

Table  36.iFffhicle  Parameters  for  Test  Case  5 
TPFiqf  dVith  Original  Time  Windows 


Type 

Aveiewe  Payload 

Daily  Cost 

K 

Beddown  Cost 

dm 

C130 

3 

10 

C-17 

9 

10 

C-5 

16 

10 

HEMTT 

y?r  7 

1 

10 

M1083 

oft  c 

1 

10 

M35 

^  8 

1 

10 

DODX 

o 

o 

60 

10 

FTTX 

ZO^DEO 

42 

10 

ITTX 

uO 

5180 

_ cm _ 

52 

10 

Table  37.  Outloading  Parameters  for  Test  Case  5 


POD 

Mode 

Outload  Capacity 

^imv 

il 

Air 

1000 

Road 

1000 

Rail 

1000 

i2 

Air 

1000 

Road 

1000 

Rail 

1000 

*Note,  for  each  POD/M 


ode  pair,  the  out: 


oad  capacity  is  assumed  constant  for  all 


days  V. 


Table  38.  Unloading  Parameters  for  Test  Case  5 


Destination 

Mode 

Unload  Capacity 

^  jmv 

Air 

1000 

Road 

1000 

Rail 

1000 

J2 

Air 

1000 

Road 

1000 

Rail 

1000 

*Note,  for  each  Destination/Mode  pair,  the  unload  capacity  is  constant  for  all  v. 
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Appendix  G.  Additional  Model  Inputs  For  Test  Case  6 


Table  39.  Outloading  Parameters  for  Test  Case  6 


POD 

Mode 

Outload  Capacity 

^imv 

il 

Air 

1000 

Road 

1000 

Rail 

1000 

i2 

Air 

1000 

Road 

1000 

Rail 

1000 

*Note,  for  each  POD/Mode  pair,  the  outload  capacity  is  assumed  constant  for  all  days 


V 


Table  40.  Unloading  Parameters  for  Test  Case  6 


Destination 

Mode 

Unload  Capacity 

^  jmv 

Air 

1000 

Road 

1000 

Rail 

1000 

J2 

Air 

1000 

Road 

1000 

Rail 

1000 

*Note,  for  each  Destination/Mode  pair,  the  unload  capacity  is  constant  for  all  v  . 
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Appendix  H.  Model  Output  Solutions  Test  Cases  2-5 


71  HEMTT{s) 

leaving  POD 

11 

for 

destination 

J1 

on 

day 

3  (ROAD) 

497.00  Short  Tons 

of  Movement 

1 

85  HEMTT{s) 

leaving  POD 

11 

for 

destination 

J2 

on 

day 

3  (ROAD) 

595.00  Short  Tons 

of  Movement 

6 

36  HEMTT(s) 

leaving  POD 

11 

for 

destination 

J1 

on 

day 

4  (ROAD) 

3.00  Short  Tons 

of  Movement 

1 

249.00  Short  Tons 

of  Movement 

2 

57  HEMTT(s) 

leaving  POD 

11 

for 

destination 

J2 

on 

day 

4  (ROAD) 

399.00  Short  Tons 

of  Movement 

7 

108  HEMTT{s) 

leaving  POD 

11 

for 

destination 

J1 

on 

day 

5  (ROAD) 

1.00  Short  Tons 

of  Movement 

2 

749.17  Short  Tons 

of  Movement 

3 

29  HEMTT{s) 

leaving  POD 

11 

for 

destination 

J2 

on 

day 

5  (ROAD) 

5.00  Short  Tons 

of  Movement 

6 

196.00  Short  Tons 

of  Movement 

8 

28  HEMTT(s) 

leaving  POD 

11 

for 

destination 

J1 

on 

day 

6  (ROAD) 

0.83  Short  Tons 

of  Movement 

3 

195.00  Short  Tons 

of  Movement 

4 

43  HEMTT(s) 

leaving  POD 

11 

for 

destination 

J2 

on 

day 

6  (ROAD) 

1.00  Short  Tons 

of  Movement 

7 

300.00  Short  Tons 

of  Movement 

9 

15  HEMTT(s) 

leaving  POD 

11 

for 

destination 

J1 

on 

day 

7  (ROAD) 

5.00  Short  Tons 

of  Movement 

4 

100.00  Short  Tons 

of  Movement 

5 

1  HEMTT(s) 

leaving  POD 

11 

for 

destination 

J2 

on 

day 

7  (ROAD) 

4.00  Short  Tons 

of  Movement 

8 

3.00  Short  Tons 

of  Movement 

10 

71  HEMTT{s) 

leaving  POD 

11 

for 

destination 

J2 

on 

day 

8  (ROAD) 

497.00  Short  Tons 

of  Movement 

10 

1  HEMTT{s) 

leaving  POD 

12 

for 

destination 

J2 

on 

day 

4  (ROAD) 

6.00  Short  Tons 

of  Movement 

14 

72  HEMTT(s) 

leaving  POD 

12 

for 

destination 

J1 

on 

day 

5  (ROAD) 

500.00  Short  Tons 

of  Movement 

11 

142  HEMTT{s) 

leaving  POD 

12 

for 

destination 

J2 

on 

day 

5  (ROAD) 

994.00  Short  Tons 

of  Movement 

14 

58  HEMTT(s) 

leaving  POD 

12 

for 

destination 

J1 

on 

day 

6  (ROAD) 

400.00  Short  Tons 

of  Movement 

12 

1  HEMTT{s) 

leaving  POD 

12 

for 

destination 

J2 

on 

day 

6  (ROAD) 

4.00  Short  Tons 

of  Movement 

15 

43  HEMTT{s) 

leaving  POD 

12 

for 

destination 

J1 

on 

day 

7  (ROAD) 

300.00  Short  Tons 

of  Movement 

13 

28  HEMTT(s) 

leaving  POD 

12 

for 

destination 

J2 

on 

day 

7  (ROAD) 

196.00  Short  Tons 

of  Movement 

15 

72  HEMTT{s) 

leaving  POD 

12 

for 

destination 

J2 

on 

day 

8  (ROAD) 

500.00  Short  Tons 

of  Movement 

16 

Figure  10.  ITDM  Solution  Test  Case  2 
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1  C130(s) 

leaving  POD 

12 

for 

destination 

J2 

on 

day 

4  (AIR) 

11.00  Short  Tons 

of  Movement 

14 

1  C130(s) 

leaving  POD 

12 

for 

destination 

J1 

on 

day 

5  (AIR) 

9.00  Short  Tons 

of  Movement 

11 

1  C130(s) 

leaving  POD 

12 

for 

destination 

J2 

on 

day 

6  (AIR) 

11.00  Short  Tons 

of  Movement 

15 

1  C130(s) 

leaving  POD 

12 

for 

destination 

J1 

on 

day 

7  (AIR) 

9.00  Short  Tons 

of  Movement 

13 

1  C130(s) 

leaving  POD 

12 

for 

destination 

J2 

on 

day 

8  (AIR) 

8.00  Short  Tons 

of  Movement 

16 

1  C130(s) 

leaving  POD 

12 

for 

destination 

J2 

on 

day 

9  (AIR) 

12.00  Short  Tons 

of  Movement 

16 

2  DODX(s) 

leaving  POD 

11 

for 

destination 

J1 

on 

day 

3  (RAIL) 

400.00  Short  Tons 

of  Movement 

1 

2  DODX(s) 

leaving  POD 

11 

for 

destination 

J2 

on 

day 

3  (RAIL) 

400.00  Short  Tons 

of  Movement 

6 

1  DODX(s) 

leaving  POD 

11 

for 

destination 

J1 

on 

day 

4  (RAIL) 

100.00  Short  Tons 

of  Movement 

1 

100.00  Short  Tons 

of  Movement 

2 

3  DODX(s) 

leaving  POD 

11 

for 

destination 

J1 

on 

day 

5  (RAIL) 

150.00  Short  Tons 

of  Movement 

2 

450.00  Short  Tons 

of  Movement 

3 

1  DODX(s) 

leaving  POD 

11 

for 

destination 

J2 

on 

day 

5  (RAIL) 

200.00  Short  Tons 

of  Movement 

6 

2  DODX(s) 

leaving  POD 

11 

for 

destination 

J1 

on 

day 

6  (RAIL) 

300.00  Short  Tons 

of  Movement 

3 

100.00  Short  Tons 

of  Movement 

4 

2  DODX(s) 

leaving  POD 

11 

for 

destination 

J2 

on 

day 

6  (RAIL) 

400.00  Short  Tons 

of  Movement 

7 

1  DODX(s) 

leaving  POD 

11 

for 

destination 

J1 

on 

day 

7  (RAIL) 

100.00  Short  Tons 

of  Movement 

4 

100.00  Short  Tons 

of  Movement 

5 

1  DODX(s) 

leaving  POD 

11 

for 

destination 

J2 

on 

day 

7  (RAIL) 

200.00  Short  Tons 

of  Movement 

8 

2  DODX(s) 

leaving  POD 

11 

for 

destination 

J2 

on 

day 

8  (RAIL) 

300.00  Short  Tons 

of  Movement 

9 

100.00  Short  Tons 

of  Movement 

10 

2  DODX(s) 

leaving  POD 

11 

for 

destination 

J2 

on 

day 

9  (RAIL) 

400.00  Short  Tons 

of  Movement 

10 

3  DODX(s) 

leaving  POD 

12 

for 

destination 

J2 

on 

day 

4  (RAIL) 

600.00  Short  Tons 

of  Movement 

14 

2  DODX(s) 

leaving  POD 

12 

for 

destination 

J1 

on 

day 

5  (RAIL) 

400.00  Short  Tons 

of  Movement 

11 

1  DODX(s) 

leaving  POD 

12 

for 

destination 

J2 

on 

day 

5  (RAIL) 

200.00  Short  Tons 

of  Movement 

14 

Figure  10.  ITDM  Solution  Test  Case  2  (cont.) 
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2  DODX(s) 

1  DODX(s) 

1  DODX(s) 
20  HEMTT(s) 
13  HEMTT(s) 
7  HEMTT(s) 
20  HEMTT(s) 
13  HEMTT(s) 
7  HEMTT(s) 
20  HEMTT(s) 
20  HEMTT(s) 


leaving  POD 

12 

for 

destination 

J1 

on 

day 

6  (RAIL) 

leaving  POD 

12 

for 

destination 

J1 

on 

day 

7  (RAIL) 

leaving  POD 

12 

for 

destination 

J2 

on 

day 

9  (RAIL) 

leaving  POD 

12 

for 

destination 

J2 

on 

day 

4  (ROAD) 

leaving  POD 

12 

for 

destination 

J1 

on 

day 

5  (ROAD) 

leaving  POD 

12 

for 

destination 

J2 

on 

day 

5  (ROAD) 

leaving  POD 

12 

for 

destination 

J2 

on 

day 

6  (ROAD) 

leaving  POD 

12 

for 

destination 

J1 

on 

day 

7  (ROAD) 

leaving  POD 

12 

for 

destination 

J2 

on 

day 

7  (ROAD) 

leaving  POD 

12 

for 

destination 

J2 

on 

day 

8  (ROAD) 

leaving  POD 

12 

for 

destination 

J2 

on 

day 

9  (ROAD) 

400.00  Short  Tons 
200.00  Short  Tons 
200.00  Short  Tons 
140.00  Short  Tons 
91.00  Short  Tons 
49.00  Short  Tons 
140.00  Short  Tons 
91.00  Short  Tons 
49.00  Short  Tons 
140.00  Short  Tons 
140.00  Short  Tons 


of  Movement 
of  Movement 
of  Movement 
of  Movement 
of  Movement 
of  Movement 
of  Movement 
of  Movement 
of  Movement 
of  Movement 
of  Movement 


12 

13 

16 

14 

11 

14 

15 

13 

15 

16 

16 


Figure  11.  ITDM  with  MPBR  Solution  Test  Case  2 
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41  C130(s) 

leaving  POD 

11 

for 

destination 

J1 

on 

day 

3  (AIR) 

492.00  Short  Tons 

of  Movement 

1 

50  C130(s) 

leaving  POD 

11 

for 

destination 

J2 

on 

day 

3  (AIR) 

600.00  Short  Tons 

of  Movement 

6 

1  C130(s) 

leaving  POD 

11 

for 

destination 

J1 

on 

day 

4  (AIR) 

8.00  Short  Tons 

of  Movement 

1 

4.00  Short  Tons 

of  Movement 

2 

33  C130(s) 

leaving  POD 

11 

for 

destination 

J2 

on 

day 

4  (AIR) 

396.00  Short  Tons 

of  Movement 

7 

83  C130(s) 

leaving  POD 

11 

for 

destination 

J1 

on 

day 

5  (AIR) 

246.00  Short  Tons 

of  Movement 

2 

750.00  Short  Tons 

of  Movement 

3 

17  C130(s) 

leaving  POD 

11 

for 

destination 

J2 

on 

day 

5  (AIR) 

4.00  Short  Tons 

of  Movement 

7 

200.00  Short  Tons 

of  Movement 

8 

16  C130(s) 

leaving  POD 

11 

for 

destination 

J1 

on 

day 

6  (AIR) 

192.00  Short  Tons 

of  Movement 

4 

25  C130(s) 

leaving  POD 

11 

for 

destination 

J2 

on 

day 

6  (AIR) 

300.00  Short  Tons 

of  Movement 

9 

9  C130(s) 

leaving  POD 

11 

for 

destination 

J1 

on 

day 

7  (AIR) 

8.00  Short  Tons 

of  Movement 

4 

100.00  Short  Tons 

of  Movement 

5 

42  C130(s) 

leaving  POD 

11 

for 

destination 

J2 

on 

day 

9  (AIR) 

500.00  Short  Tons 

of  Movement 

10 

42  C130(s) 

leaving  POD 

12 

for 

destination 

J1 

on 

day 

5  (AIR) 

500.00  Short  Tons 

of  Movement 

11 

84  C130(s) 

leaving  POD 

12 

for 

destination 

J2 

on 

day 

5  (AIR) 

1000.00  Short  Tons 

of  Movement 

14 

34  C130(s) 

leaving  POD 

12 

for 

destination 

J1 

on 

day 

6  (AIR) 

400.00  Short  Tons 

of  Movement 

12 

25  C130(s) 

leaving  POD 

12 

for 

destination 

J1 

on 

day 

7  (AIR) 

300.00  Short  Tons 

of  Movement 

13 

17  C130(s) 

leaving  POD 

12 

for 

destination 

J2 

on 

day 

7  (AIR) 

200.00  Short  Tons 

of  Movement 

15 

42  C130(s) 

leaving  POD 

12 

for 

destination 

J2 

on 

day 

9  (AIR) 

500.00  Short  Tons 

of  Movement 

16 

Figure  12.  ITDM  Solution  Test  Case  3 
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1  C5(s) 

leaving  POD 

11 

for 

destination 

J1 

on 

day 

3  (AIR) 

60.00  Short  Tons 

of  Movement 

1 

10  C5(s) 

leaving  POD 

11 

for 

destination 

J2 

on 

day 

3  (AIR) 

600.00  Short  Tons 

of  Movement 

6 

2  C130(s) 

leaving  POD 

11 

for 

destination 

J1 

on 

day 

4  (AIR) 

20.00  Short  Tons 

of  Movement 

1 

4.00  Short  Tons 

of  Movement 

2 

7  C5(s) 

leaving  POD 

11 

for 

destination 

J1 

on 

day 

4  (AIR) 

420.00  Short  Tons 

of  Movement 

1 

2  C5(s) 

leaving  POD 

11 

for 

destination 

J2 

on 

day 

4  (AIR) 

120.00  Short  Tons 

of  Movement 

7 

9  C5(s) 

leaving  POD 

11 

for 

destination 

J1 

on 

day 

5  (AIR) 

246.00  Short  Tons 

of  Movement 

2 

294.00  Short  Tons 

of  Movement 

3 

2  C5(s) 

leaving  POD 

11 

for 

destination 

J2 

on 

day 

5  (AIR) 

100.00  Short  Tons 

of  Movement 

7 

20.00  Short  Tons 

of  Movement 

8 

8  C5(s) 

leaving  POD 

11 

for 

destination 

J1 

on 

day 

6  (AIR) 

456.00  Short  Tons 

of  Movement 

3 

24.00  Short  Tons 

of  Movement 

4 

3  C5(s) 

leaving  POD 

11 

for 

destination 

J2 

on 

day 

6  (AIR) 

180.00  Short  Tons 

of  Movement 

7 

4  C5(s) 

leaving  POD 

11 

for 

destination 

J1 

on 

day 

7  (AIR) 

176.00  Short  Tons 

of  Movement 

4 

64.00  Short  Tons 

of  Movement 

5 

7  C5(s) 

leaving  POD 

11 

for 

destination 

J2 

on 

day 

7  (AIR) 

180.00  Short  Tons 

of  Movement 

8 

240.00  Short  Tons 

of  Movement 

10 

3  C130{s) 

leaving  POD 

11 

for 

destination 

J1 

on 

day 

8  (AIR) 

36.00  Short  Tons 

of  Movement 

5 

3  C17(s) 

leaving  POD 

11 

for 

destination 

J2 

on 

day 

8  (AIR) 

105.00  Short  Tons 

of  Movement 

10 

5  C5(s) 

leaving  POD 

11 

for 

destination 

J2 

on 

day 

8  (AIR) 

300.00  Short  Tons 

of  Movement 

9 

1  C17(s) 

leaving  POD 

11 

for 

destination 

J2 

on 

day 

9  (AIR) 

35.00  Short  Tons 

of  Movement 

10 

2  C5(s) 

leaving  POD 

11 

for 

destination 

J2 

on 

day 

9  (AIR) 

120.00  Short  Tons 

of  Movement 

10 

3  C17(s) 

leaving  POD 

12 

for 

destination 

J2 

on 

day 

4  (AIR) 

105.00  Short  Tons 

of  Movement 

14 

11  C5(s) 

leaving  POD 

12 

for 

destination 

J2 

on 

day 

4  (AIR) 

655.00  Short  Tons 

of  Movement 

14 

4  C17(s) 

leaving  POD 

12 

for 

destination 

J1 

on 

day 

5  (AIR) 

140.00  Short  Tons 

of  Movement 

11 

6  C5(s) 

leaving  POD 

12 

for 

destination 

J1 

on 

day 

5  (AIR) 

360.00  Short  Tons 

of  Movement 

11 

4  C5(s) 

leaving  POD 

12 

for 

destination 

J2 

on 

day 

5  (AIR) 

240.00  Short  Tons 

of  Movement 

14 

Figure  12.  ITDM  Solution  Test  Case  3  (cont.) 
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5  C130(s) 
8  C17(s) 

1  C5(s) 

5  C5(s) 

4  C17(s) 

1  C5(s) 

5  C130(s) 

4  C17(s) 

5  C5(s) 


leaving  POD 

12 

for 

destination 

J1 

on 

day 

leaving  POD 

12 

for 

destination 

J1 

on 

day 

leaving  POD 

12 

for 

destination 

J1 

on 

day 

leaving  POD 

12 

for 

destination 

J1 

on 

day 

leaving  POD 

12 

for 

destination 

J2 

on 

day 

leaving  POD 

12 

for 

destination 

J2 

on 

day 

leaving  POD 

12 

for 

destination 

J2 

on 

day 

leaving  POD 

12 

for 

destination 

J2 

on 

day 

leaving  POD 

12 

for 

destination 

J2 

on 

day 

6  (AIR) 

60.00  Short  Tons 

of  Movement 

12 

6  (AIR) 

280.00  Short  Tons 

of  Movement 

12 

6  (AIR) 

60.00  Short  Tons 

of  Movement 

12 

7  (AIR) 

300.00  Short  Tons 

of  Movement 

13 

7  (AIR) 

140.00  Short  Tons 

of  Movement 

15 

7  (AIR) 

60.00  Short  Tons 

of  Movement 

15 

9  (AIR) 

60.00  Short  Tons 

of  Movement 

16 

9  (AIR) 

140.00  Short  Tons 

of  Movement 

16 

9  (AIR) 

300.00  Short  Tons 

of  Movement 

16 

Figure  13.  ITDM  with  MPBR  Solution  Test  Case  3 
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8  C5(s) 

leaving  POD 

11 

for 

destination 

J1 

on 

day 

3  (AIR) 

480.00 

Short  Tons 

of  Movement 

1 

10  C5(s) 

leaving  POD 

11 

for 

destination 

J2 

on 

day 

3  (AIR) 

590.00 

Short  Tons 

of  Movement 

6 

4  C5(s) 

leaving  POD 

11 

for 

destination 

J1 

on 

day 

4  (AIR) 

20.00 

Short  Tons 

of  Movement 

1 

220.00 

Short  Tons 

of  Movement 

2 

2  C5(s) 

leaving  POD 

11 

for 

destination 

J2 

on 

day 

4  (AIR) 

120.00 

Short  Tons 

of  Movement 

7 

13  C5(s) 

leaving  POD 

11 

for 

destination 

J1 

on 

day 

5  (AIR) 

30.00 

Short  Tons 

of  Movement 

2 

750.00 

Short  Tons 

of  Movement 

3 

3  C5(s) 

leaving  POD 

11 

for 

destination 

J2 

on 

day 

5  (AIR) 

10.00 

Short  Tons 

of  Movement 

6 

170.00 

Short  Tons 

of  Movement 

8 

3  C5(s) 

leaving  POD 

11 

for 

destination 

J1 

on 

day 

6  (AIR) 

180.00 

Short  Tons 

of  Movement 

4 

5  C5(s) 

leaving  POD 

11 

for 

destination 

J2 

on 

day 

6  (AIR) 

280.00 

Short  Tons 

of  Movement 

7 

10.00 

Short  Tons 

of  Movement 

8 

10.00 

Short  Tons 

of  Movement 

9 

2  C5(s) 

leaving  POD 

11 

for 

destination 

J1 

on 

day 

7  (AIR) 

20.00 

Short  Tons 

of  Movement 

4 

100.00 

Short  Tons 

of  Movement 
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Figure  14.  ITDM  Solution  Test  Case  4 
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Figure  15.  ITDM  with  MPBR  Solution  Test  Case  4 
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Figure  16.  ITDM  with  MPBR  Solution  Narrow  Time  Window  Test  Case  5 
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Figure  16.  ITDM  with  MPBR  Solution  Narrow  Time  Window  Test  Case  5  (cont.) 
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Figure  17.  ITDM  with  MPBR  Solution  Wider  Time  Window  Test  Case  5 
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Figure  17.  ITDM  with  MPBR  Solution  Wider  Time  Window  Test  Case  5  (cont.) 
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Appendix  I.  MPBR  VBA  Code  Updates 


Objective  Function. 

'Build  Q  variables  (beddown  upper  bound) 

Fori  =  ITo  UBound(POD) 

For  m  =  1  To  UBound(mode) 
coeff  =  10  'g  *  3 

Print  #1, "  +  "  &  coeff  &  "  *  "  &  "Q_"  &  _ 
POD(i)  & 
mode(m); 

Next  m 
Next  i 


Print  #1, 

'Update  Status  Form 
OBJFUNC. La  belProgress. Width  =  0 

OBJFUNC.Labell. Caption  =  "Generating  Requirement  Constraints... 


Constraints. 


'*****Build  New  QjX  Constraint  ***** 

RowCount  =  1 

'initialize  current  w  value  based  off  of  RowCount=l 

currj  =  w(RowCount,  1) 

currj  =  w(RowCount,  2) 

curr_m  =  w(RowCount,  3) 

curr_k  =  w(RowCount,  4) 

daynum  =  1 
xRowCount  =  1 
LHScount  =  0 

Print  #1,  "[MinMax_"  &  constraintnum  &  "] ", 

While  RowCount  <=  UBound(w) 

If  daynum  >=  w(RowCount,  6)  And  daynum  <=  w(RowCount,  7)  Then 
While  xRowCount  <=  UBound(w) 

If  w(xRowCount,  1)  =  currj  And  _ 
w(xRowCount,  3)  =  curr_m  And  _ 

daynum  >=  w(xRowCount,  6)  And  daynum  <=  w(xRowCount,  7)  Then 
Print  #1, "  +  XJ'  &_ 
w(xRowCount,  1)  &  " J'  &  _ 
w(xRowCount,  2)  &  " J'  &  _ 
w(xRowCount,  3)  &  " J'  &  _ 
w(xRowCount,  4)  &  " J'  &  _ 
daynum; 

LHScount  =  LHScount  +  1 
End  If 

xRowCount  =  xRowCount  +  1 
Wend 
'do  RHS 

If  LHScount  >  OThen 

Print  #1,"  <=  "&"QJ'&_ 
w(RowCount,  1)  &  "  J'  &  _ 
w(RowCount,  3); 

Print  #1, 

constraintnum  =  constraintnum  +  1 
If  RowCount  =  UBound(w)  And  daynum  =  absmaxday  Then 
'nada 
Else 
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Print  #1,  "[MinMax_"  &  constraintnum  &  'print  for  next  constraint  if 
not  on  last  vehicle  constraint 
End  If 
End  If 

xRowCount  =  1 
daynum  =  daynum  +  1 
LHScount  =  0 

Elself  daynum  <  w(RowCount,  6)  Then 
daynum  =  daynum  +  1 

Else  'Implies  daynum  >  w(RowCount,  7)  Then 
daynum  =  1 

RowCount  =  RowCount  +  1 
If  RowCount  <=  UBound(w)  Then 
currj  =  w(RowCount,  1) 
currj  =  w(RowCount,  2) 
curr_m  =  w(RowCount,  3) 
curr_k  =  w(RowCount,  4) 

End  If 
End  If 

If  RowCount  Mod  10  =  0  Then 

PctDone  =  RowCount  /  UBound(w) 

With  OBJFUNC 

.FrameProgress. Caption  =  Format(PctDone,  "0%") 

.LabelProgress. Width  =  PctDone  *  (.FrameProgress. Width  - 10) 

End  With 

'The  DoEvents  statement  is  responsible  for  the  form  updating 
DoEvents 
End  If 

Wend 


'update  status  form 

OBJFUNC. LabelProgress. Width  =  0 

OBJFUNC. Labell. Caption  =  "Generating  Outloading  Constraints..." 
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MINIMAL  VEHICLE  BEDDOWN 
APPROXIMATIONS  FOR  THE  IMPROVED 
THEATER  DISTRIBUTION  MODEL 


Appendix  J.  Research  Summary  Chart 
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